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SUMMARY 


Trie  successful  design  of  a  propeller- type  v/STOL  aircraft 
requires  the  precise  calculation  and  associated  optimization 
of  the  propeller  performance  at  the  static  condition.  A  gen¬ 
eral  theory  for  performance  calculations  has  been  formulated 
based  on  a  continuous  vortex  representation  along  the  lines 
of  the  classical  lifting -line  model. 

As  opposed  to  forward  flight,  the  deformation  of  the  wake  is 
appreciable  just  behind  the  propeller,  and  its  determination 
constitutes  the  heart  of  the  static  problem. 

A  computer  program  has  been  developed  to  calculate  both  the 
inflow  at  the  propeller  and  the  induced  velocity  at  any  field 
point  for  an  arbitrary  description  of  the  trailing  vortex 
sheets.  To  approximate  the  force-free  condition  imposed  on 
the  wake,  an  initial  wake  hypothesis  derived  from  the  theory 
of  the  Generalized  Actuator  Disk  was  first  used.  The  result¬ 
ing  comparisons  with  both  detailed  and  gross  measurements 
were  unsatisfactory  and  a  refined  hypothesis  was  derived. 

The  refined  wake  hypothesis  provides  a  more  reasonable  repre¬ 
sentation  of  the  "pitch"  of  the  elements  of  the  deformed 
trailing  vortex  sheets  as  well  as  the  envelope  of  their  tra¬ 
jectories.  Using  this  hypothesis,  comparisons  at  one  power 
coefficient  show  excellent  agreement  with  test  data,  the  dis¬ 
crepancy  in  thrust  and  pitch  setting  being  about  1%  and  2$ 
respectively,  and  the  discrepancy  in  the  figure  of  merit 
being  about  . 
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INTRODUCTION 


The  design  of  high-performance  propellers  for  v/STOL  aircraft 
has  been  hindered  by  the  lack  of  a  suitable  theory  for  the 
determination  of  performance  at  the  hover  operating  condition. 
This  theory  is  needed  for  calculating  the  inflow  velocity  at 
each  blade  station  so  that  two-dimensional  airfoil  data  can 
be  applied  to  find  the  overall  thrust  and  obtain  the  design 
characteristics  of  the  blade  for  peak  overall  airplane  per¬ 
formance.  The  static  performance  of  conventional  propellers 
has  not  been  a  critical  design  condition. 

The  static  thrust  generated  by  conventional  propellers  has 
been  calculated  using  single  point  methods  based  on  full- 
scale  propeller  test  data.  These  procedures  do  not  show  how 
detailed  blade  geometry  changes  affect  the  overall  perform¬ 
ance  and,  therefore,  are  not  suitable  for  blade  design  stud¬ 
ies.  To  overcome  this  difficulty,  attempts  were  made  to 
extrapolate  the  theory  for  the  forward  flight  condition  to 
the  zero  speed  case  and  then  calculate  the  inflow  velocity. 

This  procedure  appeared  to  be  suitable  for  conventional  pro¬ 
peller  design  purposes,  but  it  did  not  have  the  necessary 
accuracy  for  use  in  designing  propellers  for  v/STOL  aircraft. 

The  need  for  peak  propeller  performance  with  V/STOL  aircraft 
is  greater  than  with  conventional  aircraft  as  the  propellers 
must  generate  a  total  lift  equal  to  the  gross  weight  of  the 
aircraft.  Since  the  payload  of  these  vehicles  is  20 %  to  30$ 
of  the  gross  weight,  a  small  improvement  in  the  lift  or  hover 
performance  of  the  propellers  results  in  a  large  improvement 
in  payload  capacity.  High  levels  of  cruise  efficiency  must 
also  be  obtained  to  keep  the  fuel  weight  requirements  to  a 
minimum . 

The  principal  difficulty  in  computing  the  inflow  velocity  in 
the  static  case  is  locating  the  trailing  vortex  elements. 

The  position  of  the  elements  in  the  wake  is  determined  by  the 
propeller  rotation  and  tdie  velocities  generated  by  the  blades; 
i.e.,  these  elements  must  deform  such  that  they  are  force  free. 
At  the  forward  flight  condition,  however,  the  wake  position 
is  determined  mainly  by  the  propeller  rotation  and  the  axial 
free  stream  velocity.  Thus,  in  forward  flight  the  position 
of  the  wake  is  known  based  on  the  operating  condition,  whereas 
in  the  static  case  the  position  of  the  wake  is  unknown.  Since 
the  inflow  is  extremely  sensitive  to  the  wake  position,  the 
deformation  of  the  trailing  vortex  elements  must  be  found  pre¬ 
cisely  for  the  elements  near  the  blade  station  at  which  the 
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inflow  is  being  evaluated.  This  is  especially  important  for 
the  outboard  stations,  where  small  changes  have  the  most 
significant  effect  on  the  performance. 

The  complicated  interaction  between  the  wake  and  the  inflow 
as  described  makes  it  necessary  to  use  the  process  of  itera¬ 
tion  to  solve  the' problem.  Thus  for  an  assumed  load  distri¬ 
bution  and  wake  deformation,  the  inflow  is  calculated  and  so 
from  two-dimensional  airfoil  data  a  new  load  distribution  is 
obtained.  With  the  new  load  distribution,  a  new  wake  defor¬ 
mation  is  computed  and  the  process  repeated  until  reasonable 
convergence  is  realized. 

Since  the  propeller  and  the  turbine  engine  are  an  effective 
lightweight  system  for  providing  both  hover  lift  and  propul¬ 
sion  at  cruise  for  V/STOL  aircraft,  the  ability  to  design 
the  optimum  configuration  is  of  considerable  importance.  To 
fill  the  existing  void,  a  theoretical  program  was  undertaken 
so  that  the  necessary  data  could  be  calculated  for  use  in 
the  design  of  propellers.  Effort  was  made  to  apply  all 
available  test  data  as  a  guide  for  the  analysis  as  well  as  a 
check  on  the  results.  The  experience  gained  from  the  devel¬ 
opment  and  use  of  the  forward  fliqht  theory  was  employed  to 
keep  the  program  practical  and  useful. 


2 


CHAPTER  1 


THEORETICAL  FORMULATION 


Basic  Equations 


We  consider  a  propeller  rotating  at  a  constant  angular  speed 
Q  in  a  uniform,  inviscid  fluid  which  has  a  density  p  and 
is  at  rest  at  infinity.  The  blade  thickness -to -chord  and 
chord-to-radius  ratios  are  assumed  to  be  sufficiently  small 
that  the  classical  lifting-line  formulation  is  a  reasonable 
approximation.  We  can  then  represent  the  propeller  blades 
by  radial  bound  vortex  lines  accompanied  by  a  force-free 
system  of  trailing  vortex  sheets. 

A  propeller-fixed  cylindrical  coordinate  system  (x,r,0)  is 
chosen  with  the  axis  of  rotation  as  the  x-axis,  see  Fig.  1  . 
The  N  blades  are  located  in  the  plane  x  =  0  and  arranged 
symmetrically  such  that 


0p  =  2tt(p-1)/N  p  «  1,  2,  ...  ,  N  (1) 


i.e.,  the  p  =  1  blade  coincides  with  0=0. 


The  induced  velocity  anywhere  in  the  propeller  field  is  de¬ 
termined  by  integration  of  the  Biot-savart  law  over  all  the 
vortex  elements  comprising  the  system  in  terms  of  their 
strength,  orientation  and  distance  to  the  field  point  of 
interest,  and  may  be  expressed  in  the  form 


U“  I  /P[r(rp^B 

P=1  R,, 

V“  I  /P[r(lp)y8 

p=l  R,, 


ar(rp) 

dr 

P 


dr(rp) 

drp 


(2) 


(3) 
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where 

u 

v 

w 

*h 

RP 

r 

-dr 

r 

P 

n 

rT 


is  the  axial  component  of  the  velocity  induced  at  a 
field  point  (x,r,9)  and  taken  positive  in  the 
positive  x -direction, 

is  the  corresponding  radial  component,  positive  in 
the  positive  r-direction,  and 

is  the  corresponding  tangential  component,  positive 
in  the  positive  6 -direction? 

is  the  radius  of  the  propeller  hub; 

is  the  radius  of  the  propeller  tip; 

is  the  bound  blade  circulation  strength  at  (0,r  ,0  ) 

Jr  Jr 

is  the  strength  of  the  vortex  sheet  element  trailing 
from  (O,rp,0p)  ; 

is  the  radial  variable  of  integration  along  the  p^ 
propeller  blade; 

is  the  influence  function  for  the  axial  component  of 
the  velocity  Induced  at  (x,r,0)  by  a  bound  radial 
vortex  element  of  unit  strength  and  unit  length  at 
(°,rp,0p)  , 

is  the  corresponding  function  for  the  radial  compo¬ 
nent,  and 

is  the  corresponding  function  for  the  tangential 
component;  and 

is  the  influence  function  for  the  axial  component  of 
the  velocity  induced  at  (x,r ,0)  by  an  arbitrarily 
deformed  vortex  sheet  element  of  unit  strength  and 
semi-inf iaite  length  trailing  from  (0,r  ,0  )  , 

Jr  JP 

is  the  cor  responding  function  for  the  radial 
component,  and 


5 


9L  is  the  corresponding  function  for  the  tangential 
T  component . 

The  influence  functions  can  be  derived  from  first  principles 
(Ref.  24,  pp.  1.10  to  1.12,  for  example )  or  obtained  from 
previous  results  (Ref.  20,  pp.  13  to  16).  They  are  summarized 
in  Table  1  where 

x  is  the  axial  coordinate  of  any  point  on  the 
v  trajectory  of  the  vortex  sheet  element  trailing 
from  (0,rp,ep)  , 

rv  is  the  corresponding  radial  coordinate,  and 

is  the  corresponding  angular  coordinate? 

uv  is  the  value  of  u  at  (xv<rv,0v)  , 

vv  is  the  corresponding  value  of  v  ,  and 

wv  is  the  corresponding  value  of  w  ?  and 

t  is  a  time  parameter  to  be  defined. 

When  a  point  on  a  blade  is  chosen  as  a  field  point,  the  con¬ 
tribution  from  the  bound  vertex  for  this  blade  is  irrelevant 
and  so  discarded,  as  in  wing  theory  (Ref.  13,  pp.  13i  to  133)- 

The  contribution  from  the  bound  vortices  to  the  induced  veloc¬ 
ity  is  simple  to  calculate.  The  contribution  to  the  induced 
velocity  from  the  elements  of  the  trailing  vortex  sheets,  on 
the  other  hand,  is  very  complicated.  Not  only  do  the  influ¬ 
ence  functions  require  an  integration  over  the  length  of  the 
element,  but  what  is  even  worse,  the  location  or  "trajectory” 
of  the  element  is  not  fixed.  Rather,  these  elements  must 
drift  force  free  by  aligning  themselves  with  the  streamlines 
of  the  flow,  which  in  differential  form  (Ref.  24,  pp.  1.1  to 
1.2,  for  example)  are 


dx  _  dr  _  d0 
u  “  v  fl+(  w/r) 


(5) 


These  equations  can  be  integrated  in  terras  of  the  parametric 
time  t  that  it  takes  a  fluid  particle  to  move  along  this 
trajectory,  giving 
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(6) 


t 

x  =  /  u  dr 

V  J  V 

0 

t 

rv  =  rp  +  f  vv  <7) 

0 

t 

ev  =  ep  +  at  +  J  (wy^)  dt  (8) 

o 

The  integrands  in  Eqs.  (6)  to  (8)  are  functions  of  the 
coordinates  along  the  trajectory  and  the  coordinates  them¬ 
selves  are  functions  of  t  . 

For  convenience,  we  now  introduce  a  set  of  operators  such 
that  Eqs.  (2)  to  (4)  become  simply 


^U(T,  -dr/dr?  u,v,w) 

(9) 

0v(r,  -dr/dr?  u,v,w) 

(10) 

0w(r,  -dr/dr;  u,v,w) 

(11) 

where 

0  is  the  Biot-Savart  operator  for  the  axial  induced 
velocity  component, 

0  is  the  corresponding  operator  for  the  radial 
component,  and 

0  is  the  corresponding  operator  for  the  tangential 
component . 

Note  carefully  that  the  operators  operate  on  the  bound  and 
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trailing  vortex  strengths  with  the  velocity  field  itself  as  a 
parameter.  These  equations  constitute  a  set  of  simultaneous, 
nonlinear,  singular  integral  equations  over  the  domain  of  the 
blades  and  the  trailing  vortex  sheets.  This  set  is  complete 
if  r  is  specified. 

If  r  is  not  specified,  another  equation  is  required  to  re¬ 
late  the  circulation  strength  and  velocity  field  through  the 
propeller  geometry.  Consider  the  force  and  velocity  diagram 
for  a  typical  section  of  the  blade  at  the  radial  station 
(0,r  ,9  )  ,  see  Fig.  2  ,  where 

-t*  Jr 

Wp  is  the  total  local  velocity  seen  by  the  blade  section; 

u  is  the  local  axial  inflow,  i.e.  the  axial  component 

P  of  the  induced  velocity  at  (0,r  ,9  )  , 

P  P 

v  is  the  corresponding  radial  inflow  (not  shown,  but 
p  perpendicular  to  the  plane  of  the  section),  and 

w  is  the  corresponding  tangential  inflow; 

r 

a  is  the  local  blade  angle  of  attack; 

£  is  the  local  blade  pitch  setting; 

c p  is  the  angle  between  ftr  and  W  ; 

p  p  p 

7  is  the  tangent  angle  of  the  sectional  lift -drag  polar; 

dD  is  the  elemental  profile  drag  on  the  blade  section;  and 

dL  is  the  elemental  lift  on  the  blade  section. 

From  the  definition  of  the  sectional  lift  coefficient  in 
terms  of  the  local  blade  chord  b  , 


CL  2  2(dL/dr  )/pW 2  b 


(12) 


the  Kutta-Joukowski  formula  (Ref.  24,  pp.  5.21  to  5-22,  for 
example)  can  be  written  as 


T  =  b  VT  C_/2 
P  >•> 


(13) 
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Figure  2.  Velocity  and  Force  Diagram  at  Typical  Section. 


where  the  lift  coefficient  is  assumed  known  (Ref.  J)  as  - 
function  of  a  for  a  given  section  at  approximately  th 
local  Mach  and  Reynolds  numbers.  This  is  the  equation 
needed  in  conjunction  with  Eqs.  (9)  to  (11)  to  complete  the 
set  when  the  propeller  geometry  is  specified. 


In  either  case,  the  circulation  or  the  propeller  geometry 
specified,  the  prediction  of  the  overall  performance  can  be 
computed  once  the  inflow  has  been  calculated.  In  particular, 
if  we  resolve  the  resultant  of  dL  and  dD  as  indicated  by 
the  dotted  lines  in  Fig.  2  and  integrate  along  the  blade  for 
the  N  blades,  the  total  thrust  T  and  power  P  in  non- 
dimensional  form  are 


where  the  thrust  and 
defined  by 


cos(0p+7)  /wp\2 

Gb  cos  7  \£2R/ 

P 

c  siptyo-)  /«P\S 

L  COS  7  VfiRp/ 

power  coefficients  1 


the  propeller  are 


p  p  2i 

CT  s  7T  T/4Pn  Rp  = 
Cp  5  AAp^R p5  = 


respectively,  with  Q  =  27m  and 


p  4 

T/pn^Dp 

(16) 

p/pn3Dp5 

(17) 

*p =  V2  • 


From  the  thrust  and  power  coefficients  we  can  calculate  the 
figure  of  merit  f/M  (Ref.  19#  pp.  353  to  355)/  given  in 
percent,  by 


F/M 


=79-8  CT3/2/cp 


(18) 
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To  conclude  this  section,  it  is  of  interest  to  discuss  the 
nondimensionalization  of  the  basic  equations  of  the  static 
thrust  problem.  As  we  have  implied  by  the  form  of  Eqs.  (16) 
and  (17)/  the  characteristic  time  is  1/Cl  and  the  character¬ 
istic  length  is  Rp  .  If  the  equations  are  nondimensionalized 
with  respect  to  these  quantities,  complete  similarity  exists 
with  the  exception  of  the  sectional  lift  and  drag  data  which 
in  general  must  include  the  effects  of  Mach  number  and 
Reynolds  number,  as  noted  earlier. 


Outline  of  solution 

The  set  of  equations  described  above  cannot  be  solved  analyt¬ 
ically,  so  approximate  numerical  techniques  must  be  employed. 
Iteration  is  perhaps  the  most  straightforward  approach  and 
the  one  which  we  have  chosen. 

First  we  rewrite  Eqs.  (9)  to  (ll)  in  a  different,  but  equiv¬ 
alent,  form,  namely 


u  =  u  -  ku[u  -  #u(r,  ...  )]  (19) 


v  =  v  -  Kv[v  -  ...  )]  (20) 


w  =  w  -  k^w  -  ^w(r,  ...  )]  (21) 


where 

K  is  a  factor  chosen  as  necessary  to  achieve  convergence 
u  of  the  axial  component  of  the  induced  velocity  field, 

K  is  the  corresponding  factor  for  the  radial  component 
v  and 

K  is  the  corresponding  factor  for  the  tangential  com- 
w  ponent . 

These  iteration  factors  may  be  constant  or  may  depend  on  any 
of  the  variables  of  the  problem. 
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With  the  propeller  geometry  specified,  iteration  begins  by 
finding  a  zeroth,  or  initial,  approximation  to  the  inflow 
and  so  the  blade  circulation  through  Eq.  (13)#  the  Kutta- 
Joukowski  formula.  The  first  approximation  to  the  inflow 
follows  by  operating  with  the  Biot-Savart  operator  on  the 
zeroth  approximation  to  the  blade  circulation  and  some  guess 
at  the  induced  velocity  field  in  the  propeller  wake.  These 
results,  along  with  the  iteration  factors  and  the  zeroth 
approximations  to  the  inflow  and  the  wake  velocity  field, 
are  then  substituted  into  the  right-hand  sides  of  Eqs.  (19) 
to  (21)  to  give  the  first  approximations  to  the  inflow  and 
the  wake  velocity  field,  respectively.  The  iteration  con¬ 
tinues  by  operating  as  before  on  the  first  approximation  to 
determine  a  second,  and  so  forth.  Once  convergence  within 
a  prescribed  accuracy  is  reached,  the  propeller  performance 
can  be  calculated. 

For  the  case  with  the  blade  circulation  specified,  on  the 
other  hand,  iteration  would  proceed  similarly  except  that  the 
circulation  would  not  change  at  each  step.  Upon  convergence, 
the  blades  could  be  designed  consistent  with  the  Kutta- 
Joukowski  formula  and  the  performance  calculated. 

The  iteration  approach  necessitates,  at  each  step  of  the 
iteration,  calculation  of  the  induced  velocity  field  at  a 
sufficient  number  of  points  to  define  accurately  the  local 
streamlines  in  the  propeller  wake.  Unfortunately,  the  number 
of  calculations  required  for  adequate  representation  of  the 
complete  distributed  vortex  sheets  is  very  large.  As  a  re¬ 
sult,  we  looked  for  some  way  which  would  reduce  the  number  of 
calculations  but  still  satisfy  the  force-free  condition  of 
the  wake,  at  least  in  an  approximate  sense. 

The  question  of  convergence  and  the  proper  choice  of  the 
iteration  factors  in  Eqs.  (19 )  to  (21)  present  a  major  dif¬ 
ficulty.  Standard  methods  to  guarantee  convergence  (Ref.  27, 
for  example)  are  net  feasible  for  calculations  of  this  scope. 
We  will  elaborate  more  on  convergence  later. 


Pitchfork  Model 


When  our  original  investigation  was  undertaken  in  October 
1962,  little  experimental  information  was  available  on  the 
detailed  nature  of  the  propeller  wake  deformation,  other  than 
a  few  smoke  pictures.  Consequently,  we  turned  to  intuitive 
considerations  for  a  plausible  wake  behavior.  Based  upon 
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order -of -magnitude  studies  using  results  from  finite  wing 
theory,  we  hypothesized  that  roll  up  of  the  trailing  vortex 
sheet  from  each  blade  into  a  pair  of  concentrated  vortices 
was  the  predominant  feature  of  the  wake  and  formulated  our  so- 
called  "Pitchfork  Model"  (Ref.  22).  This  model  treated  roll 
up  in  a  simple  fashion  by  asstiming  that  the  trailing  vortex 
sheet  from  each  blade  was  completely  undeformed  axially  and 
radially  until  the  inboard  and  outboard  roll-up  points  were 
reached,  whereupon  the  respective  portions  of  the  sheet  were 
assumed  to  contract  discontinuously,  i.e.  to  roll  up  instan¬ 
taneously,  to  a  pair  of  concentrated  vortices. 

The  Biot-Savart  operator  appropriate  to  this  model  was  pro¬ 
grammed  in  steps  for  numerical  computation  on  the  CDC  1604 
digital  computer  at  the  Cornell  University  Computing  Center 
(Refs.  22  and  8).  First  approximations  in  the  iteration 
were  calculated  at  three  pitch  settings  of  the  3(10188a2P2) 
propeller  (Refs.  9  and  10;.  At  one  of  these  settings  a  para¬ 
metric  study  of  the  effects  of  roll-up  point  variation  on  the 
first  approximation  was  made  (Ref.  9)  and  a  second  approxima¬ 
tion  was  calculated  for  the  most  promising  of  these  (Ref.  10). 
The  results  indicated  that  the  Pitchfork  Model  was  indeed  a 
step  towards  an  improved  theory,  and  gave  good  qualitative  in¬ 
dications  of  the  effect  of  overall  contraction. 

However,  quantitative  predictions  were  impaired  greatly  by 
the  failure  to  represent  the  initial  axial  and  radial  defor¬ 
mation  of  the  distributed  trailing  vortex  sheets,  especially 
near  the  tip.  We  thus  concluded  that  continuous  deformation 
of  the  vortex  sheets  was  necessary  in  an  adequate  model  of 
the  vortex  wake.  No  conclusion  about  the  importance  of  roll 
up  could  be  drawn  at  that  time,  although  the  lack  of  funda¬ 
mental  knowledge  of  the  basic  nature  of  roll  up  in  the  pro¬ 
peller  case  was  clear  as  well  as  the  need  for  more  detailed 
experimental  investigations  to  improve  the  wake  model. 


continuous  Deformation  Model 


When  the  present  theoretical  investigation  was  initiated,  we 
reviewed  extensively  previous  theories  for  propellers  in  for¬ 
ward  flight  with  the  specific  aim  of  gaining  a  better  under¬ 
standing  of  the  overall  nature  of  the  trailing  vortex  sheet 
deformation  (Ref.  6,  pp.  218  to  222;  Ref.  18;  Ref.  26,  pp.  77 
to  87;  Refs.  28  and  14).  In  general,  the  sheets  stretch 
axially,  together  with  a  radial  contraction  and  a  tangential 
distortion.  The  edges  of  the  sheets  are  also  locally  un¬ 
stable  and  tend  to  roll  up  as  for  a  finite  wing. 
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While  these  deformations  are  not  of  primary  significance  for 
the  practical  determination  of  aircraft  propeller  performance 
in  forward  flight  (Ref.  3)/  they  are  the  heart  of  the  static 
problem.  Without  a  free  stream,  not  only  do  the  deformations 
occur  much  closer  to  the  propeller  and  so  exert  more  influ¬ 
ence,  but  they  are  also  much  larger.  In  fact,  the  axial  elon¬ 
gation  which  is  only  a  perturbation  of  the  basic  helical 
sheets  in  forward  flight  is  now  the  total  variation  itself. 

Simultaneously,  we  turned  to  the  new  experimental  knowledge 
of  the  wake  that  was  available  in  the  summer  of  1964.  We  had 
taken  photographs  of  the  condensation  paths  of  the  tip  vor- 
ticity.  Also,  Canadair,  Ltd.  had  made  new  smoke  visualization 
studies . 

Based  upon  these  results  and  the  above  discussions,  we  decided 
to  concentrate  upon  the  continuous  axial  stretching  and  radial 
deformation  of  the  trailing  vortex  sheets,  or  in  other  words, 
a  "Continuous  Deformation  Model".  At  the  same  time,  we  con¬ 
cluded  that  it  was  important  to  provide  for  incorporation  of 
the  effect  of  roll  up,  if  necessary,  at  a  later  date. 

Inasmuch  as  the  associated  development  of  a  satisfactory 
approximation  to  the  force-free  condition  of  the  wake  con¬ 
tinued  throughout  the  investigation,  we  will  defer  the  details 
until  the  calculations  for  the  specific  propellers  are  de¬ 
scribed. 
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CHAPTER  2 


METHOD  OF  CALCULATION 


jlementation 


In  order  to  carry  out  the  iterative  solution  which  we  have 
outlied,  a  general  computer  program  was  developed  to  evalu¬ 
ate  the  induced  velocity  field  given  by  Eqs.  (2]  to  (4).  The 
remaining  calculations  involved,  namely  Eq.  (13)  and  Eqs. 

(19)  to  (21),  are  straightforward  and  were  carried  out  by 
hand.  Once  the  inflow  at  the  propeller  was  established,  the 
resulting  predictions  of  the  overall  performance  expressed 
by  Eqs.  (14),  (15)  and  (18)  were  found  by  means  of  a  modifi¬ 
cation  of  the  integration  method  of  the  existing  C-W  strip 
analysis  (Ref.  11).  This  also  was  carried  out  by  hand. 

In  this  chapter  we  will  describe  the  computer  program,  as 
well  as  the  associated  numerical  analyses,  and  indicate  its 
use  to  date. 


Description  of  computer  Program 


The  computer  program  was  made  as  flexible  as  possible  to  pro¬ 
vide  a  "working  tool"  for  handling  whatever  details  we  might 
wish  to  incorporate  as  the  study  progressed.  It  can  treat  an 
arbitrary,  continuous  deformation  of  the  trailing  vortex 
sheets,  with  or  without  discontinuous  roll  up  into  discrete 
vortices.  If  present,  roll  up  can  occur  either  outboard,  in¬ 
board,  or  both  outboard  and  inboard.  The  discrete  vortices 
themselves  are  assumed  to  be  regular  helices.  Their  geometry 
is  fixed  by  the  local  pitch  of  the  element  of  the  trailing 
vortex  sheet  adjoining  the  roll-up  point  but  the  location  of 
this  point  may  be  varied.  Besides  these  basic  options,  al¬ 
ternate  routes  were  included  in  the  program  to  minimize  the 
running  time  wherever  we  could. 

The  program  was  written  in  FORTRAN  63  for  the  CDC  1604  Com¬ 
puter.  The  particular  machine  that  we  used  is  located  at 
the  Cornell  University  Computing  Center.  To  achieve  the 
flexibility  desired,  most  of  the  principal  operations  were 
relegated  to  subroutines,  leaving  the  logic  for  the  assorted 
options  as  the  main  framework  of  the  program. 
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After  several  runs,  we  saw  that  it  was  important  to  examine 
the  trajectories  of  the  trailing  vortex  sheet  elements  for 
physical  reasonableness  without  running  the  full  program. 
Therefore  a  short  program  was  abstracted  to  compute  just  the 
trajectories.  It  was  generally  run  before  each  run  of  the 
full  program. 

From  the  duplicate  FORTRAN  63  decks  of  both  programs  which 
have  been  cut  and  supplied  to  USAAVLABS,  complete  listings 
can  be  made. 

In  the  following  we  will  emphasize  the  option  for  continuous 
deformation  without  roll  up,  since  this  is  the  one  which  we 
used  for  almost  all  of  our  runs.  Instructions  are  provided 
for  the  corresponding  preparation  of  input  data  in  Appendix 
I  . 


Representation  of  Circulation  Distribution 


In  general,  the  blade  circulation  distribution  r  is  tabu¬ 
lated  numerically  at  several  radial  stations,  but  for  certain 
of  the  subsequent  integrations  we  need  r  in  analytic  form. 

A  suitable  form  is  a  Glauert-type  series  (Ref.  13,  pp.  138  to 
139)/  or 


£ 

r  =  X  s±ni  ^p  (22) 

i=l 

rP  s  "  %(VV  cos  *p  (23) 


where 

G^,  is  the  i  nondimens ional  Glauert  coefficient?  and 

i>  is  the  Glauert  variable  which  runs  from  0  to  tt 

P  with  its  center  midway  between  the  hub  and  tip. 
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This  form  possesses  a  square-root  behavior  at  both  the  tip 
and  the  hub.  For  the  tip  this  is  the  proper  behavior  re¬ 
quired  for  a  lifting -line  formulation  (Ref.  5,  pp.  171  to 
174,  for  example).  For  the  hub  it  is  approximately  correct 
when  the  blade  necks  down  rapidly  at  the  shank  as  the 
3(13l68A10P3)  propeller  does.  When  the  hub  is  like  the  root 
of  a  wing,  on  the  other  hand,  the  slope  of  the  circulation 
distribution  should  be  zero  from  a  rigorous  viewpoint.  This 
could  have  been  achieved  exactly  by  an  additional  condition 
on  the  determination  of  the  Glauert  coefficients,  but  prac¬ 
tically  the  contribution  of  the  hub  to  the  overall  perform¬ 
ance  is  so  small  that  zero  slope  can  be  achieved  with  or 
without  this  condition. 

There  are  various  ways  that  the  Glauert  coefficients  can  be 
evaluated.  W<^  chose  the  method  of  least  squares  which  mini¬ 
mizes  the  error  at  the  specified  data  points  (Ref.  23/  PP* 
363  to  370).  The  program  provides  for  the  numerical  values 
of  r  to  be  given  at  up  to  30  points,  with  as  high  as 

20  .  Usually  <£  =  15  was  adequate.  For  a  fairly  smooth 
circulation,  the  worst  error  at  any  point  was  on  the  order 
of  1$,  and  for  an  irregular  shape,  about  twice  as  great. 


Velocity  Induced  by  Bound  Blade  Vortices 

As  is  generally  known,  the  contribution  of  the  bound  blade 
vortices  to  the  induced  velocity  is  identically  zero  at  the 
propeller  blades,  see  Table  1,  but  is  generally  nonzero  at 
field  points  located  off  the  blades.  Consequently,  the 
appropriate  influence  functions  are  regular  throughout  the 
range  of  interest.  This,  combined  with  the  algebraic  form 
of  these  functions  and  the  fixed  location  of  the  vortex 
elements,  makes  the  integration  of  the  first  terms  of  Eqs. 

(2)  to  (4)  relatively  simple. 

The  integration  proceeds  as  follows.  The  influence  functions 
are  calculated  at  the  radial  stations  for  which  the  data  for 
the  circulation  distribution  are  tabulated.  These  values  are 
multiplied  by  the  corresponding  values  of  r  to  give  a  set 
of  values  for  the  integrands.  Values  of  the  integrand  not 
contained  in  this  set  are  found  as  necessary  by  3-point 
Lagrangian  interpolation  (Ref.  23,  pp.  72  to  73,  for  example). 
The  integration  itself  then  is  carried  out  over  the  original 
radial  variable  by  a  generalized  Gaussian  scheme  devised  by 
A.  L.  Kaskel  at  TAR,  Inc.  In  this  scheme,  a  10-point 
Gaussian  quadrature  is  first  applied  to  the  total  interval 
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of  integration.  The  total  interval  is  in  turn  subdivided 
into  two  intervals  and  the  10-point  Gaussian  quadrature 
applied  to  each  of  these  intervals  and  the  results  added  to 
get  the  answer.  Next,  the  total  interval  is  subdivided  into 
three  intervals,  and  so  on.  At  each  step  after  the  second, 
the  answer  is  compared  with  the  two  values  which  immediately 
precede  it,  ard  the  subdivision  is  continued  until  these 
three  values  are  within  a  desired  accuracy. 

The  accuracy  used  was  +  0.00005  of  the  tip  rotational 
speed.  For  most  of  our  ruins  this  required  3  to  5  subdi¬ 
visions  of  the  total  interval,  or  equivalently,  the  weighting 
of  the  integrand  at  30  to  50  points. 


Velocity  Induced  by  Trailing  Vortex  Sheets 


As  pointed  out  earlier,  we  must  in  general  calculate  the 
trajectories  of  the  elements  of  the  trailing  vortex  sheets  in 
order  to  evaluate  the  contribution  by  these  sheets  to  the  in¬ 
duced  velocity.  Physically,  most  of  the  wake  deformation 
occurs  just  downstream  of  the  propeller,  say  within  a  dis¬ 
tance  of  a  blade  radius  or  so.  Integration  of  Sqs.  (6)  to 
(8),  therefore,  is  necessary  only  in  this  region.  Further 
downstream  the  trajectories  may  be  approximated  as  regular 
helices  with  suitable  pitch,  either  distributed  or  contracted 
discontinuously  into  discrete  vortices. 

To  carry  out  the  trajectory  integrations,  a  prescribed  veloc¬ 
ity  field  along  the  trailing  vortex  sheets,  either  from  a 
guess  or  the  preceding  step  of  the  iteration,  has  to  be 
stored  in  the  computer  as  input  data.  Though  the  real  flow 
field  is  not  axisymmetric,  we  can  in  effect  imagine  that  it 
is  since  we  only  move  along  a  trajectory.  In  other  words, 
the  flow  is  periodic  in  9  between  the  trailing  vortex 
sheets  so  that  u  ,  v  and  w  at  any  set  of  values  for  6 
which  are  27t/N  apart,  at  any  x  and  r  ,  are  identical. 
This  is  illustrated  in  Fig.  3  for  the  cross  section  of  the 
wake  of  a  3-bladed  propeller  with  the  trailing  vortex  sheets 
arbitrarily  deformed.  The  actual  radial  velocity  is  given 
by  the  solid  line  and  the  apparent  axisymmetric  radial  coun¬ 
terpart  is  given  by  the  dotted  line. 

The  stored  velocity  field  is  specified  at  a  number  of  values 
of  x  and  r  and  3-*P°int  Lagrangian  interpolation,  first  in 
r  and  then  in  x  ,  carried  out  for  other  values  in  between. 
For  each  of  15  radial  stations,  we  took  20  axial  stations 
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between  the  propeller  plane,  or  slightly  upstream  when 
necessary,  and  some  distance  downstream. 


Figure  3«  Actual  and  Axisymmetric  Counterpart 
for  Radial  Component. 


With  the  velocity  field  stored,  the  trajectory  can  be 
evaluated  in  a  step-by-step  marching  scheme  using  a  2 -term 
constant  slope  method.  The  first  step  off  the  blade  is  made 
by  multiplying  the  inflow  components  by  a  small,  but  arbi¬ 
trary,  time  increment  At  to  find  the  new  location  of  the 
vortex  sheet  element.  The  second,  and  successive,  steps  are 
taken  by  interpolating  to  determine  the  velocity  at  this  new 
location,  multiplying  it  by  2  At  and  adding  the  result  to 
the  coordinate  of  the  previous  location.  With  this  sch.me 
the  error  is  of  the  order  of  the  cube  of  At  as  opposed  to 
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the  square  of  At  ,  which  an  ordinary  2-term  Taylor  series 
would  afford.  The  increment  At  is  constant  as  the  march 
proceeds  hut  may  be  increased  to  a  new  value  at  one  point 
downstream  to  speed  the  march. 

Simultaneously  with  this  march,  it  is  convenient  to  perform 
the  integration  over  t  for  evaluation  of  this  contribution 
of  the  deformed  wake  to  the  influence  functions,  see  Table  1 
The  integrands  are  calculated  at  each  new  location  in  the 
march  and  integrated  by  Simpson ' s  rule  after  each  two  steps . 
Running  totals  of  the  influence  functions  are  carried  as  the 
march  proceeds. 

Assessment  of  this  scheme  was  made  by  calculation  of  an  ex¬ 
ample  for  which  the  trajectory  was  known  exactly.  Agreement 
between  the  exact  and  approximate  trajectories  was  excellent. 
In  addition,  a  number  of  tests  on  actual  trajectories  were 
made  with  different  values  of  At  .  It  was  found  that  the 
size  of  At  was  more  critical  for  the  integration  of  the 
influence  functions  to  a  desired  accuracy  than  for  the  tra¬ 
jectory  itself.  After  evaluation,  we  chose  I2At  =  0.008  for 
the  initial  increment  and  fiAt  =  0.104  for  the  second  incre¬ 
ment.  The  number  of  steps  at  the  former  was  fixed  at  58, 
but  the  number  at  the  latter  is  determined  by  the  time  taken 
for  the  "slowest  element"  to  reach  a  distance  downstream  as 
prescribed  by  the  wake  model,  see  Eq.  (6).  Usually  this  was 
about  100  to  200  steps. 

When  this  prescribed  distance  is  reached,  the  trajectories 
become  regular  helices  and  Eqs.  (6)  to  (8)  could  be  integra¬ 
ted  explicitly  because  the  axial  and  tangential  velocity  com¬ 
ponents  in  the  field  are  constant  and  the  radial  component  is 
zero.  This  was  not  done  in  the  present  computer  program  and 
the  contributions  of  these  helical  portions  to  the  influence 
functions  were  calculated  in  the  same  fashion  as  for  the  de¬ 
formed  portions. 

Since  the  helical  portions  extend  downstream  to  infinity,  it 
was  necessary  to  find  suitable  asymptotic  formulas  for  them. 

We  had  concluded  from  our  previous  studies  that  formal  asymp¬ 
totic  esqpansions  of  the  contributions  to  the  influence  func¬ 
tions  could  be  made  (Ref.  8,  p.  17) »  hut  an  even  better 
choice  is  the  calculation  of  the  ©-average  of  these  contri¬ 
butions  (Ref.  8,  pp.  18  to  23;  Ref.  9/  pp.  10,  19  to  20; 

Ref.  10,  p.  21).  These  contributions  to  the  influence  func¬ 
tions  can  be  extracted  from  the  results  of  the  theory  of  the 
Generalized  Actuator  Disk  (Ref.  14,  pp.  13  to  15).  Their  cal¬ 
culation  involves  simply  the  evaluation  of  some  well-defined 
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functions,  namely  the  Legendre  functions  of  the  second  hind 
and  plus  and  minus  one-half  order  and  the  Heuman  lambda 
function.  We  chose  to  evaluate  the  Legendre  functions  by 
relating  them  to  the  complete  elliptic  integrals  of  the  first 
and  second  hinds  (Ref.  4,  p.  249),  which  are  computed  using 
a  rational  approximation  (Ref.  1,  pp.  591  to  592).  For  the 
Heuman  lambda  function  we  tooh  an  exact  series  expansion 
(Ref.  4,  p.  300)  and  related  it,  also,  to  the  complete  ellip¬ 
tic  integrals  of  the  first  and  second  hind.  These  integrals 
are  calculated  as  before.  Where  applicable,  the  addition 
formula  for  the  Heuman  lambda  function  (Ref.  4,  p.  36)  is 
employed  to  speed  convergence  of  the  series. 

When  the  calculation  of  the  influence  functions  is  complete, 
we  multiply  them  by  their  corresponding  trailing  vortex 
strengths  and  proceed  to  the  integrations  of  the  second 
terms  of  Eqs.  (2)  to  (4).  Unfortunately,  the  singularities 
of  the  integrand  complicate  the  picture  and  require  special 
treatment  in  order  to  carry  out  these  integrations  quichly 
and  accurately. 

We  first  see  that  the  square-root  behavior  of  Eq.  (22)  mahes 
the  strength  of  the  elements  of  the  trailing  vortex  sheets 
infinite  at  the  hub  and  the  tip.  This  is  no  problem,  though, 
if  we  carry  out  the  integrations  in  terms  of  the  Glauert 
variable  instead  of  the  original  radial  variable,  as  in  finite 
wing  theory.  A  more  serious  difficulty  is  the  existence  of 
singularities  which  arise  in  the  influence  functions  when  the 
field  point  coincides  with  an  element  of  the  trailing  vortex 
system.  The  nature  of  these  singularities  can  be  determined 
by  direct  investigation  of  the  equations  in  Table  1  .  In 
particular,  at  the  lifting  line  itself  the  singularities  are 
not  integrable  for  a  nonzero  value  of  the  radial  inflow. 

This  is  the  same  nonintegrable  singularity  which  has  arisen 
in  attempts  to  extend  lifting-line  theory  to  swept -back  wings 
(Refs.  12,  25/  17/  for  example).  We  handled  this  diffi¬ 
culty  by  taking  the  first  two  steps  in  the  marching  scheme 
with  zero  radial  velocity  and  constant  axial  and  tangential 
velocity,  i.e.  as  regular  helical  elements,  after  which  the 
trailing  elements  deform  in  the  desired  fashion.  Since  At 
is  taken  so  small,  this  displacement  from  the  lifting  line  is 
infinitesimal  compared  to  the  blade  chord. 

As  shown  elsewhere  (Ref.  21,  pp.  14  to  1 6),  we  are  left  with 
two  types  of  singularities,  Cauchy  and  logarithmic.  To  handle 
the  integration  over  these  singularities,  we  sum  the  influence 
functions  over  all  the  blades  and  break  the  results  into 
"singular"  and  "regular"  parts.  If  we  define 
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where  u  and  w  are  the  inflow  at  r  from  the  previous 
approximation,  these  singular  parts  can  be  integrated  analyt¬ 
ically  in  the  usual  principal  value  sense  and  the  regular 
parts  can  be  integrated  by  the  generalized  Gaussian  scheme 
described  previously.  In  these  integrations,  the  derivative 
of  r  is  found  by  differentiation  of  Eq.  (22).  The  regular 
parts  are  determined  numerically  by  subtracting  the  singular 
parts  from  the  summation  over  all  the  blades  of  the  influence 
functions . 

The  computer  program  can  accommodate  up  to  108  elements  for 
the  trailing  vortex  sheet  from  any  blade.  Generally,  we  made 
our  runs  with  45  to  60  elements,  between  20  and  35  of  which 
were  taken  in  the  immediate  neighborhood  of  the  singularities 
to  obtain  adequate  definition  of  the  regular  parts.  The 
accuracy  specified  for  the  Gaussian  integration  was  again 
t  0.00005  of  the  tip  rotational  speed.  This  required  4  to 
20  subdivisions  between  the  hub  and  tip. 

If  the  induced  velocity  off  the  blades  and  trailing  vortex 
sheets  is  desired,  the  contribution  from  the  trailing  vortex 
sheets  has  no  singular  part,  and  so  integration  proceeds  ex¬ 
actly  as  for  the  regular  part. 


Use  of  Main  Computer  Program 


A  total  of  31  runs  was  made  with  the  main  computer  program, 

10  of  which  were  for  development  purposes .  These  are  summa¬ 
rized  in  Table  2  to  show  the  specific  propellers  considered 
and  pitch  settings  at  the  0.7  radius,  together  with  the  pur¬ 
pose  of  the  run.  In  the  propeller  designations,  the  digit 
in  front  (3  for  all  these  particular  propellers)  is  the  num¬ 
ber  of  blades,  while  the  digits  and  letters  within  the  paren¬ 
theses  refer  to  the  blade  design,  see  Appendix  II. 
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For  the  earlier  runs,  the  FORTRAN  63  deck  was  used.  Subse¬ 
quently,  a  BINARY  deck  was  cut  once  some  confidence  had  been 
established.  This  saved  about  6  minutes  in  compilation. 

The  running  times  for  the  BINARY  deck  varied  from  15  to  40 
minutes.  The  major  part  of  this  time  was  spent  on  the  march 
along  the  trajectories,  the  remaining  operations  taking  about 
5  minutes  in  any  run.  since  the  number  of  steps  at  the 
initial  time  increment  was  fixed,  the  variation  in  running 
time  can  be  traced  essentially  to  the  number  of  steps  required 
at  the  second  time  increment  to  reach  the  ax.ial  location  for 
termination  of  the  march.  This  depends  only  upon  the  axial 
component  of  velocity  along  the  trajectory  which  reaches  this 
termination  point  last. 

Description  of  the  wake  hypotheses  and  other  details,  as  well 
as  the  results  of  the  runs,  will  be  presented  in  the  next 
chapter . 
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CHAPTER  3 

THEORETICAL  RESULTS 


General 


As  summarized  in  Table  2,  calculations  were  carried  out  for 
four  propellers:  the  3(10l88A2P2) ,  the  3(109654),  the 
3(109652)  and  the  3(13l68A10P3) .  In  this  chapter  we  will 
present  the  specific  orientation  of  these  calculations  for 
each  of  the  propellers,  the  input  data  including  the  wake  hy¬ 
pothesis  used,  the  pertinent  results  and  the  intermediate 
conclusions  reached,  comparison  with  test  data,  where 
appropriate,  will  be  presented  and  discussed  in  Chapter  5* 


3 ( 10188A2P2 )  Propeller 

The  objectives  of  the  calculations  for  the  3(10l88A2P2)  pro¬ 
peller  were  twofold:  first,  to  check  out  the  main  computer 
program;  and  second,  to  compare  the  Pitchfork  Model  with  the 
Continuous  Deformation  Model.  Accordingly,  no  comparisons 
with  test  data  were  planned. 

The  choice  of  input  was  made  as  follows.  The  Pitchfork  Model 
described  earlier  for  finite  roll-up  times  was  rerun  at  both 
the  10.0°  and  20.0°  pitch  settings  (Ref.  9 »  p.  17*  Ref.  10,  p. 
16 ) .  Runs  were  also  made  at  these  settings  for  infinite  in¬ 
board  and  outboard  roll-up  times.  The  continuous  Deformation 
Model  was  run  only  at  the  10.0°  setting.  As  an  initial  wake 
hypothesis  or  approximation  to  the  force-free  condition  of  the 
wake,'  we  assumed  that  the  variation  of  the  Induced  velocity 
along  constant  radii,  regardless  of  the  blade  circulation, 
was  given  by  the  Generalized  Actuator  Disk  theory  (Ref.  14, 
pp.  21,  29).  This  variation  was  modified  slightly  by  assum¬ 
ing  that  the  wake  reached  its  final,  regular  helical  shape 
at  a  finite  distance  downstream,  namely  one  and  one-half  radii, 
instead  of  an  infinite  distance.  Thus,  the  wake  could  be 
calculated  completely  on  a  theoretical  basis.  The  airfoil 
data  for  these  calculations  were  found  following  the  procedure 
of  the  C-W  static  strip  analysis  (Ref.  11,  pp.  8  to  12,  19  to 
21)  with  corrections  to  a  tip  speed  of  650  feet  per  second. 

The  zeroth  approximations  to  the  axial  and  tangential  inflow 
components  also  were  based  on  this  analysis.  Since  it  does 
not  provide  the  radial  component,  we  took  for  the  Continuous 
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Deformation  Model  the  finite  roll-up  time  values  with  some 
smoothing  of  the  distribution  near  the  blade  tip.  The  itera¬ 
tion  factors  in  Eqs.  (19)  to  (21)  were  all  taken  equal  to  one. 

The  resulting  first  approximations  to  the  axial,  radial  and 
tangential  inflow  components  are  presented  in  Table  3#  to¬ 
gether  with  the  zeroth,  or  initial,  approximation.  The  Pitch- 
fork  Model  results  for  finite  roll-up  times  agreed  well  with 
the  original  results.  The  Continuous  Deformation  Model  gave 
an  axial  component  over  the  middle  of  the  blade  comparable  to 
that  of  the  Pitchfork  Model  with  finite  roll-up  times.  Both 
were  greater  than  the  zeroth  approximation  there  but  fell 
below  it  outboard  of  the  0.7  radial  station  with  the  Con¬ 
tinuous  Deformation  Model,  even  leading  to  negative  values 
near  the  tip.  On  the  other  hand,  the  infinite  roll-up  time 
values  were  considerably  greater  than  all  the  others  over  the 
entire  blade.  Investigation  of  the  contribution  of  the  ele¬ 
ments  of  the  vortex  sheet  trailing  inboard  and  outboard  of 
the  maximum  of  the  blade  circulation  showed  that  this  increase 
arose  solely  from  the  outboard  vortex  elements.  That  is,  the 
axial  velocity,  and  thus  the  pitch  of  these  elements,  was  so 
low  near  the  tip  that  the  magnitude  of  their  interference  on 
the  passing  blades  was  extremely  sensitive  to  their  precise 
length.  This  trend  with  roll-up  time  was  the  same  that  was 
seen  in  our  earlier  studies  (Ref.  9,  p.  25).  At  the  20.0° 
setting  on  the  other  hand,  the  axial  velocity  and  pitch  were 
higher  near  the  tip  and  the  finite  and  infinite  roll-up  time 
results  differed  by  less  than  5$  over  the  entire  blade  for 
all  inflow  components.  The  radial  inflow  components  at  10.0° 
generally  were  consistent  with  the  axial  components  through¬ 
out.  There  were  only  small  differences  in  the  tangential 
inflow  components,  lending  further  support  to  our  previous 
observation  (Ref.  9,  pp.  10,  20)  that  the  tangential  veloc¬ 
ity  is  a  function  only  of  r/r  (cf.  Ref.  14,  p.  17). 

These  results  showed  that  the  main  computer  program  worked 
successfully  and  that  the  envelopes  of  the  trajectories  of 
the  elements  of  the  trailing  vortex  sheets,  as  predicted 
by  the  initial  wake  hypothesis,  were  physically  reasonable. 

To  investigate  this  hypothesis  in  more  detail,  we  proceeded 
to  the  3(109654)  propeller. 


3(109654)  Propeller 

The  objectives  of  the  calculations  for  the  3(109654)  propeller 
were  also  twofold:  first,  to  iterate  to  a  solution  of  Eqs.  (13) 
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TABLE  3 


FIRST  APPROXIMATIONS  TO  INFLOW 
VARIOUS  WAKE  MODELS 
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0  -  Zeroth  Approximation 

Pf  -  Pitchfork  Model,  Finite  Roll-up  Times 

PM  -  Pitchfork  Model,  Infinite  Roll -Up  Times 

CD  -  continuous  Deformation  Model,  Initial  Wake  Hypothesis 
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TABLE  3 


'FIRST  APPROXIMATIONS  TO  INFLOW  COMPONENTS 
VARIOUS  WAKE  MODELS 
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and  (19)  to  (21)  using  the  initial  wake  hypothesis;  and 
second,  to  compare  the  resulting  performance  prediction  with 
test  data,  at  least  in  a  preliminary  fashion. 

For  input,  a  pitch  setting  of  11.3°  was  chosen  and  the  air¬ 
foil  data  were  found  as  for  the  3(10l88A2P2)  propeller,  here 
corrected  to  a  tip  speed  of  1002  feet  per  second.  The  zeroth 
approximation  to  the  axial  inflow  component  was  calculated  as 
before,  while  -Nr/47rr  was  used  for  the  tangential  inflow 
once  the  blade  circulation  had  been  determined.  The  zeroth 
approximation  to  the  radial  velocity  component  was  guessed  so 
that  the  trajectory  envelopes  were  similar  to  those  for  the 
3(10188a2P2)  propeller.  First  and  second  approximations  were 
computed  with  the  axial  and  radial  iteration  factors  equal  to 
one  and  the  successive  tangential  components  found  as  for  the 
zeroth.  Also  first,  second  and  third  approximations  were 
computed  with  all  three  iteration  factors  equal  to  one-half. 
The  third  approximation  was  computed  after  the  calculations 
for  the  3(109652)  propeller  and  involved  a  different  treat¬ 
ment  of  the  hub  geometry  that  will  be  discussed  in  the  next 
section. 

The  zeroth  approximation  to  the  axial  inflow  component,  the 
first  and  second  approximations  with  the  iteration  factors 
equal  to  one,  and  the  first,  second  and  third  approximations 
with  the  factors  equal  to  one-half  are  given  in  Table  4  . 

The  first  approximation  to  the  axial  inflow  with  the  factors 
equal  to  one  behaved  like  the  corresponding  3(10l88A2P2) 
calculation,  but  the  second  approximation  computed  from  it 
had  a  physically  unreasonable  character  which  prohibited 
continuation  of  the  iteration.  This  character  can  be  ex¬ 
plained  by  separate  examination  of  the  contributions  to  the 
axial  inflow  from  the  singular  and  regular  parts  of  the  in¬ 
fluence  functions.  Since  f>r  dominates  u  and  w  ,  the 
strength  of  the  Cauchy  singularity  in  Eq.  (24)  does  not  vary 
greatly  over  the  blade.  While  the  strength  of  the  logarith¬ 
mic  singularity  does  vary  considerably,  its  total  contribu¬ 
tion  is  much  smaller  than  that  of  the  Cauchy  singularity 
(Ref.  10,  p.  23).  Therefore  the  singular  axial  inflow  con¬ 
tribution  depends  primarily  upon  the  trailing  vortex  sheet 
strength.  In  the  first  approximation  there  was  a  heavy 
concentration  of  trailing  vorticity  outboard  of  the  0.9 
radial  station.  This  led  to  the  rapid  increase  in  the 
singular  contribution  to  the  axial  Inflow  just  inboard  of 
the  tip.  The  same  effect  was  noted  in  our  earlier  calcula¬ 
tions  with  the  Pitchfork  Model  (Ref.  10,  p«  22).  On  the 
other  hand,  the  regular  contribution  to  the  axial  inflow  is 
a  function  of  the  detailed  shape  of  the  trajectories,  as 
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SUCCESSIVE  APPROXIMATIONS  TO  AXIAL  INFLOW  COMPONENT  u/ftR 


0.0510  -0.0001  0.0712  0.0255  0.0460  -0.0009 


well  as  of  their  strength.  Here,  as  a  result  of  moving 
ahead  of  the  propeller  initially,  the  contracting  trailing 
vortex  elements  passed  very  close  to  the  next  blade  at  a 
reduced  radius;  e.g.,  the  vortex  element  trailing  from  the 
tip  passed  0.012  radii  behind  the  next  blade  at  the  0.923 
radial  station.  The  strong  interference  that  resulted  was 
observed  in  the  regular  contribution  to  the  axial  inflow. 

It,  too,  had  the  behavior  associated  with  a  concentration  of 
vorticity,  which  now,  however,  was  just  inboard  of  the  0.923 
station.  Namely,  there  was  a  maximum  in  the  regular  con¬ 
tribution  at  about  0.90,  a  zero  at  about  0.93/  and  a  minimum 
(which  is  negative)  at  about  0.95/  outboard  of  which  it  in¬ 
creased.  Summing  these  singular  and  regular  contributions 
then  gave  the  unusual  character  of  the  total  axial  velocity 
distribution.  In  order  to  provide  a  convergent  iteration  it 
was  necessary  to  consider  different  iteration  factors.  We 
chose  one-half,  recalculated  the  first  and  second  approxima¬ 
tions  accordingly,  and  obtained  a  convergent  trend.  The 
radial  inflow  in  all  approximations  varied  with  the  axial 
component  in  a  way  so  that  overall  there  were  only  small 
adjustments  in  the  trajectory  envelopes  from  step  to  step. 

The  tangential  inflow  correlated  with  r/r  as  before. 

We  gained  considerable  insight  into  the  question  of  conver¬ 
gence  of  the  iteration  from  a  simplified  theory  which  we 
derived  for  static  propellers.  This  is  an  adaptation  of  the 
Generalized  Actuator  Disk  theory  along  lines  also  put  forth 
independently  by  M.  Iwasaki  (Refs.  15  and  1 6).  Basically 
only  the  radial  variation  of  the  mean,  or  0-average,  axial 
inflow  is  considered  and  all  trailing  vortex  distortion  is 
neglected.  The  system  reduces  accordingly  to  a  single 
transcendental  equation  which  is  equivalent  to  Eq.  (9)  but 
exhibits  complete  independence  of  the  blade  sections.  To 
solve  this  equation,  we  rewrote  it  in  the  form  of  Eq.  (19) 
and  iterated  as  in  our  complete  analysis.  We  found  diver¬ 
gence  if  the  iteration  factor  were  one,  and  convergence  if 
it  were  one-half.  The  exact  conditions  for  convergence 
(Ref.  23/  PP*  184  to  187,  for  example)  confirmed  these  re¬ 
sults.  Although  our  full  system  of  equations  is  more  com¬ 
plicated  and  has  many  features  which  are  different,  still 
the  simpler  system  provided  understanding  and  guidelines  at 
least  for  choosing  the  iteration  factors  for  convergence. 

The  principal  conclusion  of  the  first  two  approximations  with 
iteration  factors  of  one-half  and  one  was  that  the  former 
value  led  towards  a  convergent  iteration,  comparison  with 
experiment  at  this  single  pitch  setting  was  inconclusive  so 
further  calculations  were  undertaken  for  the  3(109652) 
propeller. 
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3(109652)  Propeller 
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The  objectives  of  the  calculations  for  the  3(109652)  propeller 
were  threefold:  first,  to  predict  the  performance  over  a  broad 
range  cf  pitch  settings,  following  the  previous  calculations; 
second,  to  compare  the  results  with  test  data;  and  third,  to 
refine  the  wake  hypothesis  as  necessary. 

For  input,  pitch  settings  of  10.0°,  l4.C°  and  16.0°  were 
chosen  for  the  initial  wake  hypothesis  and  the  10.0°  setting 
for  the  refined  wake  hypothesis.  This  refined  hypothesis  or 
approximation  to  the  force-free  condition  of  the  wake  was 
identical  to  the  initial  hypothesis  inboard  of  the  0.5  radial 
station.  Outboard,  the  variation  of  the  axial  induced  veloc¬ 
ity  component  from  the  Generalized  Actuator  Disk  theory  was 
applied  along  streamlines  rather  than  along  constant  radii. 

The  variation  of  the  radial  induced  velocity  component  was 
then  fixed  by  assuming  th~.c  both  the  initial  and  refined  wake 
hypotheses  should  give  the  same  trajectory  envelopes.  For 
the  tangential  component  the  refined  hypothesis  was  assumed 
identical  with  the  initial  hypothesis,  since  w  is  small 
compared  to  Qx  .  The  airfoil  data  were  determined  as  before 
and  were  corrected  to  a  tip  speed  of  785  feet  per  second. 

The  zeroth  approximations  to  the  axial  and  tangential  inflow 
components  were  also  found  by  the  earlier  strip  analysis. 

Three  different  distributions  of  radial  inflow  were  assumed 
at  each  pitch  setting  and  their  trajectories  were  calculated 
by  means  of  the  short  trajectory  program.  From  these,  the 
most  promising  distribution  was  selected,  based  upon  experi¬ 
ence.  All  three  iteration  factors  were  one-half  and  the 
first  two  approximations  were  computed  for  all  three  settings. 
At  the  14.0°  setting,  third  approximations  were  calculated 
with  the  hub  geometry  treated  more  like  a  wing  tip,  as  in  the 
first  two  approximations,  and  also,  more  like  a  wing  root, 
see  Appendix  II.  Thic  latter  treatment  was  continued  not 
only  for  the  third  approximations  of  this  propeller  at  the 
10.0°  and  16.0°  settings,  but  for  the  third  approximation  of 
the  3(10965*0  propeller  at  11,3°  as  well. 

In  Table  5  are  presented  the  zeroth,  first,  second  and  third 
approximations  to  all  three  inflow  components  at  10.0°,  to¬ 
gether  with  the  corresponding  values  from  the  refined  wake 
hypothesis  which  used  the  second  approximation  as  input. 
Iteration  behaved  similarly  at  all  three  pitch  settings  and 
converged  sufficiently  well  to  the  accuracy  required  for  com¬ 
parison  of  the  performance  predictions  'ith  experiment.  The 
first  and  third  approximations  to  the  axial  inflow  distribu¬ 
tion  were  very  close,  both  in  magnitude  and  shape,  while  the 
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SUCCESSIVE  APPROXIMATIONS 
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second  approximation,  though  not  greatly  different  in  magni¬ 
tude,  was  distinctly  different  in  shape.  This  trend  was 
observed  for  the  3(109654)  propeller  too.  As  can  be  seen, 
the  axial  inflow  distribution  fell  off  near  the  tip  at  10.0°; 
however,  at  14.0°  and  16.0°  the  axial  inflow  di jtributions 
remained  high  at  the  tip,  and  both  were  of  similar  shape. 

The  radial  inflow  distributions  were  remarkably  similar  in 
shape  and  magnitude  for  all  three  pitch  settings,  while  the 
tangential  inflow  again  correlated  well  with  r/r  .  In 
Fig.  4,  the  third  approximations  to  the  inflow  components 
are  plotted  for  the  initial  and  refined  wake  hypotheses.  The 
large  differences  are  clarified  by  Fig.  5  which  shows  the 
projections  for  the  respective  trajectories  of  the  tip  ele¬ 
ment  of  one  vortex  sheet.  For  the  refined  hypothesis,  we  see 
that  the  "pitch"  is  much  lower,  and  so,  the  axial  and  radial 
inflow  much  larger.  This  is  analogous  to  the  differences  be¬ 
tween  the  inflow  in  the  Pitchfork  Model  for  infinite  and 
finite  roll-up  times,  i.e.,  a  finite  roll-up  time  provides  an 
effective  pitch  increase.  The  change  in  hub  geometry  had 
only  a  local  effect  on  the  inflow.  Outboard  of  the  0.5  sta¬ 
tion  the  difference  was  less  than  3$  • 

The  first  conclusion  reached  was  that  the  refined  wake  hy¬ 
pothesis  should  be  used  for  subsequent  calculations.  This 
was  based  primarily  on  comparisons  between  theory  and  experi¬ 
ment  which  will  be  discussed  in  Chapter  5.  Although  the 
envelopes  of  the  trajectories  are  predicted  quite  well  by  the 
initial  wake  hypothesis,  their  "pitch",  which  is  critical,  is 
described  much  better  by  the  refined  wake  hypothesis.  We  also 
concluded  that  a  different  choice  for  the  iteration  factors 
was  necessary.  Study  indicated  that  further  improvement 
could  be  gained  by  choosing  one-quarter  for  these  factors. 


3(13168a10P3)  Propeller 

The  objectives  of  the  calculations  for  the  3(13168a10P3)  pro¬ 
peller  were  twofold:  first,  to  predict  the  performance  at  one 
pitch  setting  using  the  refined  wake  hypothesis;  and  second, 
to  compare  the  results  with  test  data. 

For  input,  a  pitch  setting  of  14.0°  was  chosen.  The  airfoil 
data  were  applied  in  a  slightly  different  fashion  from  that  of 
the  earlier  static  strip  analysis.  Namely,  the  basic  section 
data  (Ref.  7)  were  read  at  the  Mach  number  and  Reynolds  num¬ 
ber  appropriate  to  the  tip  speed  of  750  feet  per  second.  Al¬ 
though  the  data  were  extended  as  before  (Ref.  11,  pp.  8,  9) 
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to  high  angles  of  attack,  the  factor  of  1.1  in  the  drag  coef¬ 
ficient  (Ref.  11,  p.  21)  was  omitted.  We  determined  the 
zeroth  approximation  to  the  axial  inflow  component  as  for  the 
other  propellers  and  used  -Nr/47rr  for  the  tangential  compo¬ 
nent  once  the  blade  circulation  was  known.  Again,  three 
different  distributions  of  radial  velocity  were  tried  and  the 
one  with  the  most  promising  trajectories  was  selected.  First, 
second  and  third  approximations  were  calculated  with  the 
axial  and  radial  iteration  factors  equal  to  one-quarter,  and 
the  successive  tangential  components  were  found  as  for  the 
zeroth . 

In  Table  6  are  presented  the  zeroth,  first,  second  and  third 
approximations  to  all  three  inflow  components.  The  axial  in¬ 
flow  changed  smoothly  from  step  to  step.  Inboard  near  the 
hub  and  outboard  of  the  0.8  radial  station  the  iteration  con¬ 
verged.  Over  the  middle  of  the  blade,  the  convergence  trend 
was  inconclusive  as  far  as  the  iteration  was  carried  out. 

This  behavior  was  distinctly  different  from  our  previous  iter¬ 
ations.  The  radial  inflow  component  became  steadily  more 
negative,  which  led  to  more  realistic  trajectory  envelopes 
inboard.  As  found  throughout,  good  correlation  with  r/r  was 
observed  for  the  tangential  inflow. 

The  main  conclusion  which  we  drew  from  these  calculations  was 
that  the  convergence  properties  of  the  iteration  with  the  re¬ 
fined  wake  hypothesis  and  iteration  factors  of  one-quarter 
are  different.  It  appears  that  one  or  possibly  two  further 
approximations  are  required  to  achieve  complete  convergence 
over  the  whole  blade. 
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TABLE  6 


SUCCESSIVE  approximations  TO  INFLC 
REFINED  WAKE  HYPOTHESI 


3(13168a10P3)  Propelle 
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0  -  Zeroth  Approximation 
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2  -  Second  Approximation 

3  -  Third  Approximation 


TABLE  6 


E  APPROXIMATIONS  TO  INFLOW  COMPONENTS 
REFINED  WAKE  HYPOTHESIS 


3(13l68A10P3)  Propeller 
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CHAPTER  4 


TEST  MEASUREMENTS 


General 


It  is  desirable  to  describe  some  of  the  experimental  investi¬ 
gations  made  in  support  of  the  analysis.  These  test  data 
proved  useful  not  only  to  guide  the  theory  toward  its  final 
form,  but  to  assess  the  calculated  predictions. 

The  available  test  data  fall  into  two  categories.  The  first 
is  propeller  wake  velocity  distributions  along  a  blade  radius 
measured  with  a  hot-film  anemometer.  The  second  category  is 
overall  propeller  performance  presented  in  nondimens ional 
form  as  thrust  coefficient,  pitch  setting  and  figure  of  merit, 
all  as  a  function  of  the  power  coefficient. 

The  three  propellers  that  were  tested  are  the  3(109654),  the 
3(109652)  and  the  3(13l68A10P3) .  Blade  characteristics  for 
these  propellers  are  shown  in  Appendix  II  .  To  date,  wake 
surveys  have  only  been  made  behind  the  3(109654)  and 
3(109652)  propellers. 

All  of  the  measured  data  for  the  3(109654),  3(109652)  and 
3(13l68A10P3)  propellers  incorporated  in  this  report  were 
obtained  on  the  Curtiss  Static  Test  Rig,  which  was  placed  in 
operation  3  years  ago. 


Description  of  Curtiss  Static  Test  Rig 


The  rig  is  mounted  out  of  doors  on  a  concrete  slab,  see 
Fig.  6  .  The  propeller  is  cantilevered  well  away  from  the 
R-3350  engine  mounting  to  minimize  wake  blockage  and  inter¬ 
ference  effects.  It  is  capable  of  testing  propellers  up  to 
18  feet  in  diameter,  supplying  as  much  as  2500  horsepower. 

A  few  feet  upstream  of  the  engine  output  shaft,  the  shaft  is 
broken  to  include  a  Baldwin-Lima-Hamilton  A-l60  Torque  Cell 
as  shown  in  the  schematic  of  Fig.  7  •  Proceeding  along  the 
shaft  toward  the  propeller,  the  shaft  is  again  broken  to  con  ■ 
tain  a  metal  diaphragm.  The  purpose  of  the  diaphragm  is  to 
preload  the  tapered  roller  thrust  bearings  so  as  to  eliminate 
unwanted  deflections  of  the  bearings  as  well  as  any  thermal 
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Figure  6.  Curtiss  Static  Test  Rig. 


growth  of  the  shaft.  This  Isolates  the  thrust  load  at  the 
strain  gaged,  four -leg  thrust  ring  which  is  located  just  aft 
of  the  propeller.  The  gage  installation  is  designed  to  be 
insensitive  to  changes  in  temperature  caused  by  frictional 
energy . 

Thrust,  torque,  propeller  RPM  and  wind  velocity  signals  are 
read  out  on  four  channels  of  a  DYMEC  Digital  Data  Acquisition 
System.  Thrust  and  torque  signals  taken  from  the  thrust  ring 
and  torque  cell  are  scanned  and  printed  on  tape.  Readout  ac¬ 
curacy  is  on  the  order  of  0.1$  .  The  propeller  RPM  is  meas¬ 
ured  with  a  magnetic  pickup  directly  on  the  propeller  shaft 
and  is  also  printed  out.  Because  of  the  effect  of  wind  on 
the  performance  of  a  static  propeller,  wind  velocity  meas¬ 
ured  with  a  cup  anemometer  is  also  monitored.  Useful  data 
are  not  obtainable  when  the  steady  velocity  is  above  3  miles 
per  hour.  The  total  scanning  time  for  a  cycle  of  all  four 
channels  is  a  little  over  1  second. 

The  DYMEC  system's  output  print-out  is  in  volts,  making  it 
necessary  to  calibrate  the  system.  The  thrust  is  calibrated 
immediately  after  each  run  with  the  use  of  a  deadweight.  The 
deadweight  of  6100  pounds  is  passed  over  a  pulley  and  hooked 
to  the  propeller  hub.  A  hydraulic  ram  raises  the  weight  off 
the  ground.  The  entire  calibration  takes  about  2  mimites. 

The  torque  cell  is  calibrated  periodically  by  the  use  of 
levers  and  arms  and  requires  half  a  day  for  completion. 

Repeat  runs  on  various  propellers  have  shown  excellent  agree¬ 
ment,  and  the  data  are  felt  to  have  better  than  +  2$  accuracy. 


Description  of  Hot-Film  Anemometer 


It  is  difficult  to  obtain  accurate  propeller  wake  surveys 
since  both  the  magnitude  and  the  direction  of  the  local  ve¬ 
locity  vector  must  be  found.  In  addition,  the  instrumentation 
must  be  capable  of  measuring  the  rapid  velocity  variations 
with  time  (as  seen  here  in  space-fixed  coordinates  in  contrast 
to  the  steady  flow  seen  in  propeller-fixed  coordinates),  that 
are  associated  with  the  finite  blade  number.  This  requires  a 
very  short  response  time. 

The  instrument  used  is  the  Lintronic  Model  40W  Linear  con¬ 
stant  Temperature  Hot-Film -Anemometer.  It  is  similar  in 
principle  to  the  hot-wire  anemometer.  That  is,  a  thir.  ^re- 
heated  platinum  element,  which  is  extremely  sensitive  to 
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temperature  differentials  caused  by  variations  in  airflow,  is 
mounted  on  a  Pyrex  glass  wedge.  The  hot-film  anemometer  is 
rugged  <  :id  maintains  stability  in  physical  characteristics 
for  extended  periods  of  time,  yielding  good  signal  to  noise 
ratio  and  very  high  frequency  response.  Most  important,  it 
has  been  designed  to  yield  a  linear  variation  of  velocity 
with  electrical  signal  so  that  a  quick  calibration  at  one 
point  is  sufficient.  A  complete  description  including  cir¬ 
cuit  diagrams  can  be  found  in  the  operating  manual  (Ref.  2). 


Figure  8.  Hot-Film  Anemometer  Probe. 


The  test  probe  is  mounted  on  a  boom  that  can  be  positioned 
at  any  distance  aft  of  the  disk  and  at  any  station  along  the 
blade  radius.  The  boom  is  shown  in  the  photograph  of  Fig.  8  . 
By  remote  control  the  probe  can  be  rotated  about  its  tip  and 
also  about  its  axis  of  symmetry.  These  degrees  of  freedom 
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enable  the  element  of  the  probe  to  be  positioned  perpendicu- 
lar  to  the  resultant  velocity  vector,  much  in  the  same  way  a 
total  pressure  tap  opening  must  be  aligned  to  the  flow  to  get 
a  meaningful  result. 

The  test  probe  is  calibrated  by  dropping  the  probe  into  a 
small  test  section  upstream  of  a  blower.  This  test  section 
is  positioned  so  that  as  a  propeller  wake  is  being  surveyed, 
the  probe  can  be  immediately  dropped  into  the  test  section 
for  a  check  calibration  in  a  matter  of  seconds.  This  small 
wind  chamber -blower  arrangement  is  shown  at  the  bottom  of  the 
photograph  of  Fig.  8  .  In  the  wind  chamber,  uniform  flows  up 
to  220  feet  per  second  can  be  maintained  at  the  test  section. 
By  throttling  the  flow,  while  noting  probe  voltage,  a  cali¬ 
bration  of  voltage  versus  velocity  is  obtained. 

The  test  procedure  is  as  follows.  The  orientation  of  the 
probe  is  adjusted  until  the  maximum  signal  voltage  is  ob¬ 
tained.  Then  the  space  coordinates,  in  this  case  three 
angles,  are  recorded  together  with  the  probe  voltage  which  is 
highly  damped  to  read  an  average  value.  An  oscillograph  re¬ 
cords  the  unsteady  portion  of  the  velocity  field.  This  pro¬ 
cedure  of  maximizing  the  average  velocity  component  simpli¬ 
fies  considerably  the  test  procedure.  It  is  not  necessary  to 
seek  a  maximum  unsteady  velocity,  since  the  unsteady  portion 
of  the  velocity  is  less  than  10jfc  of  the  resultant.  Limited 
tests  have  shown  that  the  average-to-»peak-velocity  ratio  does 
not  vary  with  small  excursions  of  the  probe  alignment. 

A  curve  showing  average  axial  velocity  measurements  for  the 
3(109652)  propeller  is  presented  in  Fig.  9  •  The  probe  was 
positioned  10  inches  behind  the  trailing  edge  of  the  pro¬ 
peller  blade,  giving  a  minimum  safe  clearance  of  a  few  inches 
between  the  propeller  blade  and  the  probe  support  structure. 
An  example  of  the  data  taken  from  the  oscillograph  records  is 
shown  in  Fig.  10  .  A verage-to-peak -velocity  ratio  is  plotted 
versus  blade  radius.  These  data  were  used  to  convert  the 
data  of  Fig.  9  to  compare  with  the  inflow  at  the  blade  pre¬ 
dicted  by  the  theory.  Repeatability  of  valid  velocity  data 
is  within  £’  in  the  axial  velocity  direction. 
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CHAPTER  5 

COMPARISON  OP  THEORY  WITH  TEST  DATA 


Scope 

As  indicated  earlier,  the  theoretical  calculations  for  the 
3(10l88A2P2)  propeller  were  for  check-out  and  preliminary 
evaluation  purposes  only  and  will  not  he  compared  with  test 
data  at  all.  Since  calculations  for  the  3(10965*0  propeller 
were  also  primarily  for  theoretical  purposes,  their  compari¬ 
son  will  he  limited  in  extent.  The  principal  comparisons 
which  we  will  make  are  for  the  3(109652/  and  3(13l68A10P3) 
propellers . 


3(109654)  Propeller 

The  third  approximation  using  the  initial  wake  hypothesis 
with  the  iteration  factors  equal  to  one-half  is  chosen  as 
the  final  theoretical  result  for  the  3(109654)  propeller  at 
the  11.3°  pitch  setting. 

Comparison  of  the  axial  inflow  over  the  middle  portion  of 
the  hlade  where  data  are  available  indicates  that  the  theo¬ 
retical  values  are  on  the  order  of  65$  of  the  measured  values. 
If  the  first  or  second  approximations  are  used  instead,  these 
results  are  varied  by  about  10$  . 

The  performance  predictions  fall  roughly  7$  below  the  meas¬ 
ured  power  coefficient,  3$  above  the  measured  thrust  coeffi¬ 
cient  and  10$  above  the  measured  figure  of  merit.  Again  the 
maximum  variations  in  these  performance  predictions  in  the 
iteration  are  1$,  4$,  and  3$»  respectively. 


3(109652)  Propeller 

The  third  approximation  using  the  initial  wake  hypothesis 
with  the  iteration  factors  equal  to  one-half  is  chosan  as  the 
final  theoretical  result  for  the  3(109652)  propeller  at  the 
10.0°,  14.0°  and  16.0°  pitch  settings.  For  the  preliminary 
evaluation  of  rhe  refined  wake  hypothesis,  the  single  approx¬ 
imation  at  the  10.0°  setting  is  considered. 
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Comparison  of  the  theoretical  and  measured  distributions 
of  the  axial  inflow  at  the  14.0°  pitch  setting  is  presented 
in  Pig.  11  .  The  measured  values  were  extrapolated  to 
the  propeller  blades  from  the  data  measurements  10  inches 


-  TEST  CURVE,  CURTISS  STATIC  TEST  RIG 

■  THEORETICAL  POINT,  INITIAL  WARE  HYPOTHESIS 


behind  the  propeller,  see  Figs.  9  and  10  .  The  initial  wake 
hypothesis  was  used,  since  at  this  distance  downstream  it  does 
not  differ  much  from  the  refined  wake  hypothesis.  As  for  the 
3(109654)  propeller  the  theoretical  results  fall  generally 
below  the  measured  values  but  not  by  as  much.  At  the  0.7 
station,  the  corresponding  error  in  a  is  about  20#  .  Again, 
the  maximum  variation  of  ax^al  inflow  in  the  iteration  is 
approximately  10#  . 
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The  performance  predictions  over  the  range  of  pitch  settings 
are  given  in  Pigs.  12  to  14  .  The  overall  comparison  based 
upon  the  initial  wake  hypothesis  is  unfavorable,  generally 
the  predictions  being  worst  at  the  10.0°  pitch  setting.  Most 
disturbing,  the  entire  qualitative  shape  of  the  figure  of 
merit  versus  power  coefficient  curve  is  incorrect.  Maximum 
variations  in  the  iteration  of  the  power  and  thrust  coeffi¬ 
cients  and  figure  of  merit  are  about  2$,  4#  and  6#,  respec¬ 
tively,  at  10.0°;  2$,  2$  and  3#  at  14.0°;  and  3#/  3#  and  1# 

at  16.0°  .  The  effect  at  14.0°  of  the  different  treatment 
of  the  hub  geometry  was  much  less  than  1$  . 

The  performance  predictions  of  the  refined  wake  hypothesis 
are  also  presented  in  Pigs.  12  to  14  for  iteration  factors 
of  both  one-quarter  and  one-half.  With  either,  the  trend  is 
very  encouraging.  Furthermore,  study  indicates  that  the  re¬ 
finement  would  also  improve  the  performance  prediction  at  the 
14.0°  and  16.0°  settings.  These  results  based  on  this  single 
approximation  should  not  be  construed  as  a  final  answer  but 
are  subject  to  iteration. 


The  third  approximation  using  the  refined  wake  hypothesis  and 
iteration  factors  equal  to  one-quarter  is  chosen  as  the  final 
theoretical  result  for  the  3(l3l68A10P3)  propeller  at  the 
14.0°  pitch  setting. 


Comparisons  of  the  performance  predictions  with  the  test 
data  in  Figs.  15  to  17  are  excellent.  With  regard  to  con¬ 
vergence,  the  predictions  generally  improved  with  each  step 
in  the  iteration,  moving  monotonically  from  the  zeroth  ap¬ 
proximation  to  the  third  approximation.  The  respective 
values  found  for  the  thrust  coefficient  from  the  zeroth, 
first,  second  and  third  approximations  were  0.158,  0.151# 

0.145  and  0.140  .  For  the  power  coefficient,  the  corres¬ 
ponding  values  were  O.O638,  0.0612,  O.O587  and  0.0576  ;  and 
for  the  figure  of  merit,  they  were  78.5##  76.5#/  75*0#  and 
72.6#  .  Consequently,  we  again  see  that  further  approximations 
are  required,  as  noted  earlier. 


Summc 


of  Performance  comparisons 


Table  7  gives  an  overall  picture  of  the  performance  compari¬ 
sons  for  the  3(109654),  3(109652)  and  3(13168a10P3)  propellers, 
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Comparison . 
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comparison. 
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PITCH  SETTING  AT  O.TRn  RADIUS 


3(13168aiop3)  propeller,  flRp  =  750  fps 
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Comparison o 


Figure  17 .  Figure  of  Merit 
Comparison . 
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The  theoretical  values  are  for  the  3(10965*0  propeller,  the 
third  approximation  using  the  initial  wake  hypothesis  and 
iteration  factors  equal  to  one-half;  for  the  3(109652)  pro¬ 
peller,  the  first  approximation  using  the  refined  wake 
hypothesis  and  iteration  factors  equal  to  one-quarter;  and 
for  the  3(13168A10P3)  propeller,  the  third  approximation 
using  the  refined  wake  hypothesis  and  iteration  factors 
equal  to  one-quarter. 

As  pointed  out  earlier,  the  comparisons  for  the  3(10965*0 
propeller  are  unsatisfactory  but  those  for  the  3(109652) 
propeller  are  quite  encouraging,  although  tentative.  While 
the  comparisons  for  the  3(13168a10P3)  propeller  are  ex¬ 
cellent,  further  approximations  are  still  required. 
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CONCLUSIONS 


Based  on  the  results  obtained  in  this  study,  we  draw  the 
following  conclusions: 

A  model  based  on  continuous  deformation  of  the  wake  is 
necessary  to  calculate  the  inflow,  and  so  the  performance,  of 
a  v/sTOL  propeller  at  the  static  condition.  The  envelopes 
of  the  trajectories  of  the  elements  of  the  trailing  vortex 
sheets  and  the  "pitch"  of  these  elements  must  be  predicted 
accurately.  An  iteration  procedure  is  used  to  solve  the 
problem. 

An  initial  wake  hypothesis,  or  approximation  to  the  force- 
free  condition  of  the  wake,  gives  a  reasonable  prediction  of 
the  trajectory  envelopes,  but  an  inadequate  prediction  of  the 
"pitch".  However,  a  refined  wake  hypothesis  appears  to  give 
an  adequate  prediction  of  both  the  envelopes  and  the  "pitch"-. 

While  the  performance  predictions  for  the  3 ( 13168A10P3 )  pro¬ 
peller  based  on  the  refined  wake  hypothesis  are  excellent, 
further  approximation s)  in  the  iteration  procedure,  as  well 
as  predictions  at  other  pitch  settings  and  for  other  pro¬ 
pellers,  are  required.  Comparison  should  also  be  made  be¬ 
tween  the  wake  velocities  used  in  the  refined  wake  hypothesis 
and  those  found  both  by  actual  calculations  and  test  measure¬ 
ments  . 

The  main  computer  program  which  has  been  developed  is  very 
general.  It  can  be  utilized  over  a  range  of  forward  speeds 
and  for  a  variety  of  applications  other  than  the  present 
problem.  One  such  application  which  has  been  made  i3  the 
calculation  of  input  information  for  propeller-induced 
vibration  of  ship  hulls. 
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RECOMMENDATIONS 


Based  on  the  results  obtained  in  this  study,  we  make  the 
following  recommendations: 

Develop  additional  subroutines  to  incorporate  the  existing 
propeller  airfoil  data  into  the  main  computer  program  and 
introduce  simplifications  to  shorten  the  time  of  preparation 
and  running. 

Make  calculations  of  the  induced  velocity  in  the  wake  to 
examine  in  detail  the  refined  wake  hypothesis,  or  approxima¬ 
tion  to  the  force-»free  condition.  Make  corresponding  full 
scale  flow  measurements  to  confirm  the  theoretical  develop¬ 
ment.  When  this  hypothesis  is  established  or  corrected, 
calculate  solutions  for  three  propellers  at  three  pitch 
settings  and  compare  with  test  data. 

Calculate  the  inflow  for  four  representative  blade  load 
distributions  to  provide  an  improved  means  for  preliminary 
blade  design  estimates. 

Determine  the  optimum  circulation  distribution  for  the 
static  case  analogoue  to  the  classical  forward  flight 
optimum. 


60 


BIBLIOGRAPHY 


1.  Abzamowitz,  M. ,  and  stegun,  I.  A.,  Editors,  Handbook  of 
Mathematical  Functions  With  Formulas,  Graphs  and  Mathe¬ 
matical'  Tables,'  National  Bureau  of  standards,  Applied 
Mathematics  Series  55 ,  Washington,  D.  c.,  1964. 

2.  Anon.,  Operating  Manual  for  Linear  constant -Temperature 
Hot-Film  Anemometer  Mode£  40W,  Lintronic  Laboratories, 
Ithaca,  N.  Y.,  1963. 

3.  Borst,  H.  V.,  and  Glover,  L.  s.,  Application  of  Theodor- 
sen's  Theory  to  strip  Analysis  Procedure  for  single  "Rota¬ 
tion  Propellers,  Report  No.  C-2070,  Curtlss-wrlgnt  Corp., 
Propeller  Div.,  Caldwell,  N.  J.,  September  2,  1949- 

4.  Byrd,  P.  F.,  and  Friedman,  M.  D.,  Handbook  of  Elliptic 
Integrals  for  Engineers  and  physicTits,  sprlnger-VefliTg , 
Berlin,  Germany,  1954. 

5.  Durand,  w.  F.,  Editor,  Aerodynpiie  Theory.*  Vol.  II,  Dover 
Publications,  Inc.,  New  YGrK,  New  York,  1963. 

6.  ibid.,  vol.  IV. 

7*  Enos,  L.  H.,  and  Borst,  H.  V.,  Propeller  Performance 

Analysis  Aerodynamic  Characteristics',  NACA  l'6  Series  Air- 
lolls.  Report  No.  C-2000,  curtiss-wright  corp..  Pro- 
peller  Div.,  Caldwell,  N.  J.,  Deceittoer  2,  1948. 

8.  Erickson,  J.  c.,  Jr.,  and  Ordway,  D.  E.,  A  New  Approach  to 
the  Static  Thrust  Problem,  Phase  II,  Calculation  of  Thrust 
ancProrgue,  Model  X-lOQ  Propeller,  Sn  =  l0.0°  ,  THERM, 
Incorporated,  TAR.-TR  632,  May  1963 .  u*  'Kp 

9.  Erickson,  J.  C.,  Jr.,  and  Ordway,  D.  E.,  A  New  Approach 
to  the  static  Thrust  Problem,  Phase  III,  Further  Inflow 
Calculations  with  Extensions  to  Partially  Negative  r 
Distributions,  THERM,  Incorporated,  TAR-TR  635,  August 

19 *r. 

10.  Erickson,  J.  c.,  Jr.,  and  Ordway,  D.  E.,  A  New  Approach 
to  the  Static  Thrust  problem.  Phase  IV,  Further  Devel¬ 
opment  of  Performance  Calculation’  Method,  THERM,  Incor¬ 
porated,  TAR-TR  638,  December  1^63 • 


61 


i 


11. 


12. 


13* 


sf 


l 


£ 


14. 


15. 

16. 

17. 

18. 

19. 

20. 

21. 

22. 


Foss,  R.  L . ,  Strip  Analysis  Method  for  Evaluating 
Static  Thrust,  Report  No.  C-2402,  Curtis  s-Wrlglht 2orp . , 
Propeller  Div. ,  Caldwell,  N.  J.,  November  19,  1952. 

Gail,  A.,  Letter  to  the  Editor,  Journal  of  the  Aeronau¬ 
tical  Sciences,  volume  10,  Number  8',  October  1943, 
p.  320. 

G'.auert,  H.,  The  Elements  of  Aerofoil  and  Airscrew 
Theory,  Second  Edition,  Cambridge  University  press, 
Cambridge,  England,  1948. 

Ik  .jh,  G.  R.,  and  Ordway,  D.  E.,  The  Generalized  Actua¬ 
tor  Disk,  Therm  Advanced  Research,  Inc.,  TAR-TR  6401, 
Janvary  1964;  also  "The  Generalized  Actuator  Disk", 
Developments  in  Theoretical  and  Applied  Mechanics, 
Volume  II,  Proceedings  of  the  Second  Southeastern  con¬ 
ference  on  Theoretical  and  Applied  Mechanics,  Pergamon 
Press,  c.rford,  England,  1965,  PP.  317-336. 

Iwasaki,  M.,  Calculation  of  Static  Thrust  of  propel¬ 
lers  (1)  ,  private  Communication,  October  18,  1964. 

Iwasaki,  M.,  Calculation  of  static  Thrust:  of  a  Propel¬ 
ler  (II ) ,  Private  Communication,  December  26,  1964. 

Lin,  C.  C.,  Letter  to  the  Editor,  journal  of  tiie  Aero¬ 
nautical  sciences,  volume  11,  Number  $,  July  1944, 
pp.  195-196. 

Lock,  C.  N.  H.,  The  Application  of  Goldstein's  Theory 
to  the  Practical  Design  of  Airscrews,  British  ARC, 

R&M  No.  1377,  November  1930. 

von  Mises,  R.,  Theory  of  Flight,  Dover  Publications, 
Inc.,  New  York,  New  York,  1959. 

Ordway,  D.  E.,  Sluyter,  M.  M.,  and  sonnerup,  B.  U.  0., 
Three-Dimensional  Theory  of  Ducted  Propellers,  THERM, 
Incorporated,  TAR-TR  602,  August  I960. 

Ordway,  D.  E.,  and  Greenberg,  M.  D.,  General  Harmonic 
Solutions  for  the  Ducted  Propeller,  THERM,  Incorpora- 
ted,  TAR-TR  613/  August  ”1961. 

Ordway,  D.  E.,  and  Erickson,  J.  c.,  Jr.,  A  New  Ap¬ 
proach  to  the  Static  Thrust  Problem,  Phase  I,  Formula¬ 
tion  and  Preliminary  Assessment  of  Mathematical  Model, 
THERM,  Incorporated,  TAR-TR  6^1,  March  1963* 


62 


23.  Scarborough,  J.  B.,  Numerical  Mathematical  Analysis, 
The  Johns  Hopkins  Press,  Baltimore',  Maryland,'  1930'. 


-T*  1 

f 


24. 


25-  Theilheimer,  F.,  Letter  to  the  Editor,  journal  of  the 
Aeronautical  Sciences,  Volume  10,  Number  8 ,  October 
19537  ppT"  320-321. 

26.  Theodor sen,  T.,  Theory  of  Propellers,  McGraw-Hill,  New 
York,  New  York,  1948. 

27.  Traub,  J.  F.,  Iterative  Methods  for  the  Solution  of 
Equations ,  Prentice-Hall,  Inc.,  Englewood  cliffs, New 
Jersey,  1964. 

28.  Wu,  T.  Y.,  "Flow  Through  a  Heavily  Loaded  Actuator  Disc", 
schiffstechnik.  Band  9 >  Heft  47 ,  1962,  pp.  134-138. 


Sears,  W.  R. ,  Theoretical  Aerodynamics,  Part  1:  Intro 
duction  to  Theoretical  Hydrodynamics,  1957  Edition, 
Ithaca,  New  York,  1957- 


63 


APPENDIX  I 

OPERATING  INSTRUCTIONS  FOR  COMPUTER  PROGRAMS 


General 


Operating  instructions  for  both  the  main  computer  program 
( CODEFVEL )  as  used  for  the  continuous  Deformation  Model  and 
the  short  trajectory  computer  program  (VORTTRAJ)  are  pre¬ 
sented.  These  consist  of  descriptions  of  the  input  data  re¬ 
quired,  their  punch  format  and  the  print-out.  It  is  assumed 
that  the  basic  input  data  of  the  run  are  Known:  namely,  the 
distributions  of  the  blade  circulation,  the  inflow  velocity 
components  and  the  velocity  components  along  the  trailing 
vortex  sheets.  All  angles  must  be  in  radians  and  the  remain¬ 
der  of  the  data  must  be  nondimensionalized  with  respect  to 
the  characteristic  time  l/Q,  and  the  characteristic  length 


Main  Program  (CODEFVEL) 

Specify  the  number  of  propeller  blades  (NB),  the  propeller 
hub  radius  (RBH),  the  total  number  of  field  points  (NTIP) 
desired,  and  the  axial  (XB),  radial  (RB)  and  azimuthal  (TB) 
coordinates  of  the  field  points.  The  field  points  cannot 
exceed  20  and  must  be  numbered  consecutively  for  identifi¬ 
cation,  starting  from  that  point  on  a  propeller  blade  nearest 
the  hub  and  proceeding  along  the  blade,  followed  by  points 
off  the  blades. 

Use  the  following  values:  for  (IDIS),  3;  for  (t),  (TIR), 

(TOR)  and  (DTREG ) ,  all  0.0;  for  (INDEX),  1;  anu  for  (NREG) 
and  (IDVP),  both  0  . 

Specify  for  the  blade  circulation  distribution  the  number  of 
Glauert  coefficients  (nc),  the  maximum  value  of  the  circula¬ 
tion  (GD)  and  the  values  of  the  circulation  (gb)  at  the  total 
number  (NG)  of  radial  stations  (RG).  NC  cannot  exceed  20  and 
NG  cannot  exceed  30  . 

Specify  the  absolute  value  of  accuracy  (TOL)  required  in  the 
Gaussian  integration  scheme  and  the  maximum  number  of  stab- 
divisions  (NLIMIT)  permitted  before  termination. 

Noamnalize  the  velocity  components  at  points  along  the  trailing 
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vortex  sheets  with  respect  to  the  corresponding  values  in  the 
propeller  plane  at  the  same  radius.  Specify  the  normalized 
axial  (UB),  radial  (VB)  and  tangential  (wb)  velocity  compo¬ 
nents  at  15  radial  stations  (r)  for  each  of  20  ->,:ial  stations 
(x).  Specify  the  magnitudes  of  the  axial  (UD),  radial  (vd) 
and  tangential  (WD)  inflow  components  at  the  propeller  plane 
for  the  same  15  radial  stations. 

Specify  the  step  sizes  (DTOIC  and  DTICDH)  of  the  first  and 
second  time  increments  in  the  marching  scheme,  the  number  of 
steps  at  each  (NIOIC  and  NIICDH)  and  the  ratio  (NR)  of  DTICDH 
to  DTOIC.  NR  must  be  an  integer.  Determine  the  time  elapsed 
during  the  march  (TDH)  by  adding  the  product  of  NIOIC  and 
DTOIC  to  the  product  of  NIICDH  and  DTICDH.  Set  the  quan¬ 
tities  TDHI  and  TDHO,  the  elapsed  time  (TA)  after  which  the 
asymptotic  formulas  of  the  influence  functions  are  picked  up, 
and  the  quantities  TAI  and  TAO  all  equal  to  the  value  of  TDH. 
Specify  the  absolute  value  of  accuracy  (TOLS)  required  in 
computing  the  Heuman  lambda  function  (use  0.0001). 

Specify  the  total  number  of  trailing  vortex  elements  (NVP) 
required  for  computing  the  integrands  of  Eqs.  (2)  to  (4)  and 
their  initial  radial  coordinates  (RBP)  as  they  leave  the 
blades.  Set  their  initial  axial  (XBP)  and  azimuthal  (TBP) 
coordinates  equal  to  0.0  .  Number  the  vortex  elements  con¬ 
secutively  from  1,  which  must  coincide  with  the  hub,  to  NVp, 
which  must  coincide  with  the  tip.  NVP  cannot  exceed  108  . 
Specify,  for  each  vortex  element,  the  total  number  of  field 
points  (NIP)  at  which  the  influence  functions  for  this  ele¬ 
ment  are  desired  and  the  index  numbers  of  these  points 
( INIP ) .  specify  correspondingly,  for  each  field  point,  the 
total  number  of  vortex  elements  (NIVP)  contributing  influence 
functions  and  the  index  numbers  of  these  elements  (IVP). 

Check  NIP,  INIP,  NIVP  and  IVP  for  consistency. 

Specify  the  total  number  of  vortex  elements  (NPVP)  at  which 
intermediate  trajectory  and  influence  function  print -out  is 
desired  and  the  index  numbers  of  these  elements  (IPVP).  If 
none  is  desired,  set  NPVP  equal  to  1  and  IPVP(l)  equal  to  0 

The  punch  format  for  program  CODEFVEL  is  given  in  Table  8  and 
includes  the  input  quantities,  the  field  specification  and 
the  number  of  cards  required. 

The  print -out  of  program  CODEFVEL  is  comprehensive  and  con¬ 
sists,  in  general,  of  the  following:  the  axial,  radial  and 
tangential  velocity  components  in  the  stored  input  field? 
identification  of  the  field  points?  the  Glauert  coefficients 
and  values  for  the  representation  of  the  circulation  and  its 
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TABLE  8 


INPUT  PUNCH  FORMAT  -  PROGRAM  CODEFVEL 


Field 

No.  Of 

Input  Quantity 

Spec. 

X(l) 

X(10)  \ 

F8 

X{11) 

X(20)  / 

2 

R(l) 

R(10)  \ 

F8 

R(X1) 

R(15)  f 

2 

UB  (XI,  Rl) 

-  UB  (XIO,  Rl) 

UB (Xll,  Rl) 

-  U3  (X20,  Rl) 

. 

'  F8 

30 

UB  (XI,  R15) 

-  UB  (XIO, R15) 

UB  (Xll,  R15) 

-  UB (X20, R15) 

VB  (XI,  Rl) 

-  VB  (XIO,  Rl) 

VB  (Xll,  Rl) 

-  VB  (X20,R1) 

• 

1  F8 

30 

VB  (Xl/R15) 

-  VB (X10/R15) 

VB  (Xll,  R15) 

-  VB  (X20,  R15) 

WB(X1,R1) 

-  WB  (XIO,  Rl) 

WB  (Xll,  Rl) 

-  WB  (X20,  Rl) 

*  I 

• 

F8 

30 

WB  (XI,  R15) 

-  WB  (XIO,  R15)  ' 

WB  (Xll,  R15) 

-  WB (X20, R15)  , 

NC,  NTIP,  NVP,  NB,  NLIMIT, 

NG,  IDIS  15 

1 

UD  (1) 

UD  (10)  \ 

UD (11) 

UD (15)  j 

VD  (1) 

VD(10)  ( 

F8 

VD(ll) 

VD (15)  [ 

b 

WD(1) 

WD (10)  1 

- 

WD  (11) 

WD (15)  1 
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TABLE  8  -  contd 


TDHI,  TDHO,  TAI,  TAO,  GD, 

RBH,  TOL,  \ 

F8 

1 

TOLS 

,  TIR,  TOR  / 

<L 

XB(1) 

XB  (10)  \ 

F8 

* 

XB (11) 

XB(NTIP)  / 

RB  (1) 

RB (10)  \ 

F8 

★ 

RB  (11) 

RB(NTJP)  / 

TB  (1) 

TB  (10)  \ 

F8 

* 

TB (11) 

TB(NTIP)  ) 

GB  (1) 

GB  (10)  \ 

F8 

* 

GB (11) 

GB(NG)  / 

RG(1) 

RG(10)  \ 

F8 

* 

RG (11) 

RG(NG)  J 

NPVP,  IPVP(l) 

-  IPVP (15)  1 

L 

15 

* 

IPVP (16)  - 

IPVP  (NPVP)  J 

r 

DTOIC,  DTICDB 

DTREG,  T, 

TDH,  TA 

F8 

1 

INDEX,  NIOIC, 

NIICDH,  NR, 

NREG,  IDVP 

15 

1 

XBP,  RBP,  TBP 

- 

F8 

NIP,  INIP(i) 

-  INIP (15) 

1 

15 

** 

INIP  (16)  - 

INIP (NIP) 

/J 

NIVP,  IVP(l) 

-  IVP (15) 

l 

15 

*** 

IVP (16)  -  IVP(NIVP) 

/ 

*  -  As  Many  as  Required 

**  -  One  Set  of  as  Many  as  Required  for  Each  Vortex 

Element 

***  -  One  Set  of  as  Many  as  Required  for  Each  Field  Point 
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derivative  with  respect  to  the  Glauert  variable?  the  strengths 
of  the  Cauchy  and  logarithmic  singularities  as  defined  in 
Eqs.  (24)  to  (26)?  the  successive  values,  as  the  Gaussian  in¬ 
tegration  scheme  subdivides,  of  the  contribution  of  the  bound 
blade  vortices  to  the  induced  velocity  components  at  each 
field  point  off  the  blades;  the  contribution  from  the  bound 
blade  vortices  to  the  induced  velocity  components?  the  regular 
part  of  the  integrand  in  Eqs.  (2)  to  (4)  from  the  trailing 
vortex  sheets?  the  successive  values,  as  the  Gaussian  integra¬ 
tion  scheme  subdivides,  of  the  contribution  of  the  regular 
part  to  the  induced  velocity  components  at  each  field  point? 
the  sum  of  the  contributions  from  the  regular  part  and  from 
the  bound  blade  vortices  to  the  induced  velocity  components; 
the  contribution  from  the  Cauchy  singularity  of  the  influence 
functions  of  the  trailing  vortex  sheets  to  the  induced  veloc¬ 
ity  components?  the  contribution  from  the  logarithmic  singu¬ 
larity  to  the  induced  velocity  components?  the  resulting  total 
singular  contribution  to  the  induced  velocity  components?  and 
the  final  values  for  the  induced  velocity  components.  The  in¬ 
termediate  print-out  for  each  specified  vortex  element  con¬ 
sists  of  the  following:  the  coordinates  and  velocity  compo¬ 
nents  at  every  other  step  of  the  march  along  the  trajectory; 
the  same  quantities  at  the  termination  of  the  march?  the  con¬ 
tribution  to  the  influence  functions  at  the  specified  field 
points  by  integration  to  the  end  of  the  march?  the  total  in¬ 
fluence  f mictions  at  the  specified  field  points?  and  the 
regular  parts  of  these  influence  functions. 


Trajectory  program  (VORTTRAJ) 

specify  (UB),  (VB),  (WB),  (X),  (R),  (UD),  (VO).  (WD),  (DTOIC), 
(DTICDH),  (NIOIC),  (NIXCDH),  (nr),  (NVP),  (xbp),  (RBP)  and 
(TBP)  in  the  same  fashion  as  for  program  CODEFVEL  .  specify 
further  the  total  number  of  sets  (NVF)  of  inflow  fields  UD, 

VD  and  WD  that  are  desired.  NVF  cannot  exceed  20  and  here 
NVP  cannot  exceed  100  . 

The  punch  format  for  program  VORTTRAJ  is  given  in  Table  9  and 
includes  the  input  quantities,  the  field  specification  and 
the  number  of  cards  required. 

The  print-out  of  program  VORTTRAJ  consists  of  the  following: 
the  axial,  radial  and  tangential  velocity  components  in  the 
stored  input  field?  and  at  every  other  step  of  the  march 
along  each  trajectory,  the  coordinates,  including  the  pro¬ 
jection  (r  /R  )  cos  /  and  the  velocity  components. 

v  xr  v 


TO 
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TABLE  9 


-  PROGRAM  VOR'TTRAJ 


Input  Quantity 

Field 

Spec. 

No.  Of 
Cards 

NVP,  NVF,  NIOIC,  NIICDH,  NR  15 

1 

X(l) 

X(ll) 

X(10)  \ 

X(20)  f 

F8 

2 

R(l) 

R(ll) 

R(10)  \ 

R(X5)  / 

F8 

2 

UB  (XI,  Rl) 

UB  (Xllf  Rl) 

-  UB  (X10,  Rl)  1 

-  UB  (X20,  Rl)  i 

:  1 

\  F8 

30 

UB  (XI,  R15) 
UB  (Xll,  R15) 

-  UB  (X10,  R15)  1 

-  UB  (X20,  R15)  J 

VB  (XI,  Rl) 

VB  (Xll,  R15) 

-  VB  (X10,  Rl)  1 

-  VB (X20, Rl) 

• 

»  F8 

30 

VB  (XI,  R15) 
VB  (Xll,  R15) 

-  VB  (X10,  R15) 

-  VB  (X20,  R15)  , 

WB  (XX,  Rl) 
WB(X11,R1) 

-  WB  (XIO,  Rl) 

-  WB  (X20, Rl) 

• 

*  F8 

30 

WB(Xl,R15) 
WB  (Xll,  R15) 

-  WB (XIO, R15) 

-  WB (X20, R15) 

TABLE  9  -  contd. 


UD(R1,1)  -  UD (RIO, 1)  \ 

UD(R11,1)  -  UD (R15, 1)  / 

VD (Rl, 1)  -  VD(R10, 1) 

VD(R11, 1)  -  VD (R15, 1)  [ 

WD (Rl,  1)  -  WD  (RIO,  1) 

WD(R11,1)  -  WD(R15,1)  I 

0 

F8  * 

UD  (Rl,  NVF)  -  UD  (RIO,  NVF)  \ 

UD (Rll, NVF)  -  UD (R15, NVF)  / 

VD  (Rl, NVF)  -  VD (RIO,  NVF)  I 
VD (Rll, NVF)  -  VD(R15,NVF)  [ 

WD (Rl, NVF)  -  WD (RIO, NVF) 

WD (Rll, NVF)  -  WD (R15, NVF)  j 

XBP(l),  RBP(l),  TBP(l) 

•  F8  NVP 

XBP(NVP),  RBP(NVP),  TBP(NVP) 

DTOIC,  DTICDH  F8  1 


APPENDIX  II 


BLADE  CHARACTERISTICS 


All  three  propellers  considered  in  the  text  for  both  theo¬ 
retical  calculations  and  comparison  with  test  data,  namely 
the  3(109654) ,  3(109652)  and  3(13168a3 0P3)  propellers,  have 
three  blades,  as  indicated  by  the  digit  3  in  front  of  the 
designations.  The  digits  and  letters  within  the  parentheses 
refer  to  the  blade  design.  The  geometries  and  characteris¬ 
tics  of  these  three  blade  designs  are  given  in  Figs.  18,  19 
and  20,  with  h/b  as  the  local  blade  thickness-to-chord 
ratio  and 


where 

AF  is  the  activity  factor  of  the  blade, 

IC  is  the  integrated  design  lift  coefficient  of  the 

^i  blade  and 

C_  is  the  blade  sectional  design  lift  coefficient 
^i  (Ref.  7/  PP*  2  to  5/  for  example). 
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